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NUMERICAL METHOD AND FORTRAN PROGRAM FOR SOLVING POISSON'S 
EQUATION OVER AXISYMMETRIC REGIONS OF A SPHERE 
by Oliver W. Reese 
Lewis Research Center 

SUMMARY 

This report describes a numerical method and FORTRAN program for solving 
Poisson’s equation in axisymmetric regions of a sphere. The associated matrix equa- 
tion is solved by a direct (noniterative) method. Two geometries are considered: 

(1) A full sphere with boundary values specified on the surface only. 

(2) A spherical cone with additional boundary values specified on the conical surface 
and, optionally, along the interior negative polar axis. 

In developing the numerical method, a curvilinear cell is used at each mesh point to 
derive the finite difference equation. The points are then ordered so as to produce a 
block tridiagonal symmetric matrix when the difference equation is written at all points. 
Finally, the particular properties of this matrix are used to write an algorithm from 
which the computer program is produced. 

The program will accept various combinations of boundary conditions and spherical 
shapes. Irregular mesh spacing is permitted in both the r and 9 directions. The 
program is very easy to use and general enough to be used as a f ’black box” where the 
solution to Poisson’s equation is needed in spherical regions. The numerical method 
is also general enough for easy extension to other types of equations and other geom- 
etries. 


INTRODUCTION 

In designing a multielectrode spherical collector for a microwave amplifier tube, 
part of the problem required the solution to Poisson's equation for the electric potential 
distribution in the collector. This report describes the numerical method and the asso- 
ciated FORTRAN computer program for solving this equation in axisymmetric regions 
of a sphere. Two geometries are considered: 



(1) A full sphere with boundary values prescribed on the surface only. 

(2) A spherical cone, where additional boundary values are prescribed on the coni- 
cal surface and, optionally, along the interior negative polar axis. 

These two shapes are illustrated in figure 1. 

The following approach was used in developing the numerical method. First, a 
curvilinear cell was used at each mesh point to derive a finite difference analog to 
Poisson’s equation. Next, the mesh points were ordered in such a manner as to produce 




(a) Full sphere. 


(b) Spherical cone, 

0 O <0 <7r. 


Figure 1. - Spherical geometries considered. 


a block tridiagonal symmetric matrix when the finite difference equation was written at 
all points in the region. Finally, the particular properties of this type of matrix were 
used advantageously to produce a fast, efficient, and compact computer program. A 
direct rather than iterative method is particularly desirable in this problem since the 
program is used iteratively as a subroutine by a much larger calling program. This 
calling program continually calculates and feeds the source term used by the subroutine. 

The program is flexible in two ways: 

(1) It will accept various combinations of boundary conditions and spherical shapes. 
This feature was needed to determine the optimum design for an efficient collector. 

(2) Irregular mesh spacing is permitted in both the r and 9 directions. This 
flexibility is necessary since the solution of the problem described was nonuniform in 
character. There were regions where it varied slowly and regions where it varied 
rapidly. The area of particular interest in this study was immediately adjacent to and 
along the negative polar axis where the electron beam was entering the collector. Here, 
a relatively small mesh spacing was required to obtain the solution in sufficient detail. 
Since the use of such a small mesh throughout the entire region would be wasteful, both 
of computer storage and execution time, an irregular mesh of increasing coarseness in 
both directions was used in moving away from this axial region. Such gradation was 
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also desirable in maintaining a reasonably smooth transition from mesh point to mesh 
point in the solution. 

In addition to the features just described, the program is easy to use and general 
enough to recommend it for use as a black box in any situation where the solution to 
Poisson's equation is needed in spherical regions. 

The numerical method is also general enough for easy extension to other geometries 
and equations where a fast, compact program is needed. 

The report first explains the numerical method by discussing the derivation of the 
finite difference equation, the mesh point ordering, the generation of the tridiagonal 
block matrix, and the development of the algorithm used. 

Finally, a general description of the FORTRAN computer program is given together 
with instructions on its use. The appendixes show such details as listing of the FOR- 
TRAN program, the results of sample problem execution, the flow charts, the FORTRAN 
symbol list. 


SYMBOLS 


(All quantities are dimensionless unless otherwise specified. ) 


A i’ B i> c i 

dS 

dV 

F 

f 

G i’ K i 

L, R, T, B, D 

M 

n 

r 

Ar, A 0, A <p 
r 


r l’ r 4 


tridiagonal elements of M 
element of area 
element of volume 

coefficient of source term in Poisson's equation 
source term in Poisson's equation 
column vectors 

left, right, top, bottom, and central coefficients of the unknown at any 
general mesh point 

block diagonal matrix 

outward unit normal 

radius vector in spherical coordinates 

incremental change in spherical coordinates 

average value of r 

radius vector at midpoint of faces 1 and 4 of curvilinear cube 
unit radius vector in spherical coordinates 
total surface of curvilinear cube 
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-©-» -e- 


S r ,S 0 ,S 0 faces of curvilinear cube pierced by r, 9, and 0 axes, respectively 

AS r , ASg, AS^ areas of faces of curvilinear cube pierced by the r, 9, and 0 axes, 
respectively 

AS^, , AS„ areas of faces 1 and 4 of curvilinear cube 

r l r 4 

AS fl , AS n areas of faces 2 and 5 of curvilinear cube 

e 2 0 5 


u 

W u 0 

u 


rr 




electric potential, V 

partial derivatives of u with respect to r, 9, and 0 
a 2 u/ar 2 

3 2 u/30 2 

gradient of u taken at midpoint of faces 1 and 4 of curvilinear cube 
gradient of u taken at midpoint of faces 2 and 5 of curvilinear cube 
volume 


V i’ W i’ Z i 

e 


1 

P 

9 


e 

h 

d 2’ 0 5 


Subscripts: 

i 

j 

l 


matrices 

dielectric constant, F/M 
azimuthal angle in polar coordinates 
unit 0 vector in spherical coordinates 

O 

volume charge density, C/m 
polar angle in spherical coordinates 
average value of 9 

unit 9 vector in spherical coordinates 

average values of 6 at midpoint curvilinear cube faces 2 and 5 
vector differential operator 
Laplacian operator 

i^ ring 

j n (9 = constant) ray 


m 

n 


4 


upper limit of i 
upper limit of j 

order of block tridiagonal matrix 



STATEMENT OF PROBLEM 


The electric potential u satisfies Poisson's equation 


V 2 U = - £ 


e 


( 1 ) 


where p is the space charge density and e is the dielectric constant or permitivity. 
In spherical coordinates equation (1) is written 


or 


r 2 3r 



1 


sin 8 


_ 3 _ 

de 



1 9 2 u _ _ p 

r 2 sin 2 0 3$ 2 e 


( 2 ) 


a 2 u + 2 3u + _i_ a 2 u 
3r 2 r Sr r 2 d8 2 


cot 6 du 
r 2 38 


l a 2 u 

r 2 sin 2 0 3<p 2 



(2a) 


Azimuthal symmetry is assumed in this problem (du/dcp = 0); therefore, the fifth 
term vanishes. The particular regions of interest are bounded by the coordinates r < 1 


r\ 

, \ (a) Case I; full sphere (0q = 0). Boundary 

u 

^ Spherical surface 


Conical 

surface 



(b) Case II ; spherical segment (0q t 0). 
Boundary values are specified on spherical 
surface, and on conical surface (0q 
- constant). 


Negative 

polar 

axis-^ 



(c) Case III; spherical segment (0 q£ 0) with 
axial boundary values. Boundary values are 
specified on spherical surface, on conical 
surface (0q = constant), and on negative 
polar axis (0 = 7r line). 


Figure 2. - Three spherical configurations for which solution can be ob- 
tained. Right half of r - 0 plane shown only. Boundary conditions are 
indicated by solid lines. 
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and 0 ^ 0q < d ^ n. These regions are shown in figure 2 and are described as follows: 

Case I (0q = 0): The volume is a full sphere and boundary conditions are prescribed 
on the surface only. 

Case II (0 q + 0): The volume is a spherical segment, and additional boundary values 
are specified on the conical surface (0q = constant). 

Case HI: This case is identical to case n, with the additional option of specifying 
boundary values along the negative polar axis (0 = jt). Maximum allowable r - 0 mesh 
size is 15 by 25. Mesh size and spacing is chosen by the user. 

The mesh for the full sphere configuration was composed of (1) one center point, 

15 interior r = constant rings, and the outer or surface boundary ring; (2) twenty-five 
0 = constant lines. 


METHOD OF SOLUTION 

Of the various methods known for solving this type problem, the one presented here 
requires the least computing time and storage. 

These advantages are a consequence of ordering the mesh points in such a way as to 
produce a block tridiagonal matrix. The particular properties of this matrix are then 
used to develop a direct (noniterative) solution. 


Derivation of Finite Difference Equation 

We will now derive the finite difference equation for a general interior point by 
applying Gauss’ divergence theorem to Poisson’s equation in the cell. The resulting 
derivatives are then approximated by using a linear combination of each mesh point and 
its four neighbors. 

Poisson’s equation 


V 2 u =f 


may be written for any volume V such that 



V 2 u dV 


/ 

* / t r 


f dV 


(3) 


(4) 
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By the divergence theorem 



Vu dV = 



n 


• Vu dS 


( 5 ) 


where n is the outward unit normal and S is the surface of V. Equation (4) may now 
be rewritten as 



Vu dS = f f dV 


(6) 


Now consider a typical mesh point in spherical coordinates (r, 9, <p) enclosed by an 



Figure 3. - Curvilinear element of volume in spherical coordi- 
nates. 


element of volume as shown in figure 3. The surface of this cell is made up of three 
parts S r + Sg + S^, where S r denotes left and right faces 1 and 4, denotes top and 
bottom faces 2 and 5, denotes front and rear faces 3 and 6. The area formulas of 
these faces and the volume of the element are then 
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( 7 ) 


AS r = r 2 sin 9 AO A</> 




AS 0 = r sin 9 At A (p 

AS , = r Ar A0 

0 

dV = (r sin 9 A0)(r A0)(Ar)| 


= r sin 9 Ar A 9 A0 


Equation (6) may now be approximated by the volume bounded by this cell in the 
following manner: 


Vu = F,u r + - + 

r 




( 8 ) 


where Fp Fp and 0^ are unit vectors in spherical coordinates. Therefore, 


f- u r on face 1 


+ u r on face 4 


u fl 

— on face 2 
r 


, u 

n • Vu =-<+ — — on face 5 
r 


u 


0 


7 


r sin 9 


on face 3 


u 


0 


on face 6 


r sin 9 


(9) 


The 0 terms in expressions (8) and (9) will be neglected since u^ = 0 in this problem. 


8 



The following terms are now defined: 


2 

AS„ = r, sin 9 A9 A<b = area of face 1 
r i 1 


AS = ri sin 9 A9 Ad> = area of face 4 
r 4 4 ^ 


AS 0 = rsin # 2 Ar A0 = area of face 2 
2 


ASg = r sin Ar A<p = area of face 5 
5 


( 10 ) 


Substituting equations (8) and (9) into equation (6), we obtain 


u AS - u AS + - u„ AS. - - u . AS. = f ■ dV 
r 4 r 4 r l r l r d 5 r d 2 °2 


(U) 


where 


u 


3u 
' 1 3r 


u 


■ 

4 3r 


3u 

u n = — 

3^ 


3u 

ii. = — 

30 


r=r. 


r=r , 


0 —Be 


6~6c 


y 


( 12 ) 


Figure 4 shows that the area approximations in equation (11) may be improved by 
using the actual dimensions of the cell. 
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01 = 0 



Figure 4. - Typical portions of mesh arrangement in r - 0 
plane. Heavy lines represent front face of curvilinear 
element of volume enclosing mesh points. (1<K l, 
Kj^m). 


-2 


A V 


r i _ 1 (sin 0. + sin 0._ 1 )(0j - 0j_ 1 )A <p 


AS. 


r?(sin 0. + sin 0j_ 1 )(0j - 0j_ 1 )A0 


AS, 


(r. . sin 0. . + r. sin 0. .) _ 

= -L-k Li I Li_ ( r . - r i _ 1 )A0 


(13a) 


AS, 


(r- i sin Si + r. sin 0.) _ _ 

- — _ — z 1 — i— (r. - v. 


( r i - r i _ 1 )A <p 


J 

where 


1 < i < Z 
1 ^ j =s m 
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Finally, by substituting equations (13a), (13b), and (13c) in equation (11) and divid- 
ing through by A <p, we obtain the finite difference equation for the typical interior mesh 
point. 
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'Vi,j ~ u i,n - r* 


Ar i 


r i< e i 


? ]-i> 


r. sin 0. + r. sin 0- ■« 
J *■ J ~ _ 


where 


and where 


*i-U ~ U i,3\- 


+ l ~ "" A r i 


'r.^singj + r^sin 




sin 0j + r. j sin 0^ 


. 4 — U* *1 ( 1 • 0X11 t/: 1 T X * 0X11 L/» -I 

iiizi -v. ^ Li ill ill J - 1 


r. sin 0. i + r. ■, sin 6. 


+ [ ~ =H(r A - r A _ j )| 


r i a9 j-i 


■ < r i - >W 


r i< 9 ] - ^-i) + r,-i(8) - <>i_ j ) 


X 


fr^ sin 0^ + r i sin 0j_j + r^j sin 0^ + r i _ 1 sin ^j_A 


(14a) 


2 si < l-l 
2 < j s m 

f = - P (14b) 

£ 


Reformulation of Finite Difference Equation 

To simplify these expressions, we now use the following notation. Let 
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_ _ _ (r.i sin 0. + r.. sin 9. j) 

L i,i - r i-i (g i • e i-i> — — 

(r, sin 0, + r. sin 9 , ■. ) 


T i,j = < r i - r i-i> 


(rj sin + r i _ 1 sin 0j_j) 

( 2r i A0 i-l) 


_ _ (r* sin 9; + r. , sin 0 4 ) 

E i,i -< 7 i- 7 i-i> — L 


y as) 




r i (e i - e i-i> + r i-i (e i ~ Vi* 


r A sin 6>j + r i sin sin 0j + sin 


where 


2 < i < Z-l 


2<js m 


Substituting these expressions in equation (14a) then yields 

L i,j (u i-l,j " u i,j ) + R i,j (u i+l,j " u i, j ) + T i,j (u i,j-1 ' u i,j ) + B i,j (u i,j+1 ' u i,j ) = F i,j (f) 

(16) 


or 


L i,j u i-l,j + R i, j u i+l, j + T i,j u i,j-1 + B i,j u i,j+1 " u i,j (L i,j + R i,j + T i,j + B i,j ) f ’ F M 

(17) 


13 



Finally, 


D. .u. . + L. -u. 


+ R, .u. 


+ T. ,u. 


+ B, -u. 


i, H, j T i> j i~l> 3 ■" i> 3 i+l> j T A i, T "i,H, j+1 


= F. -f 

1,1 


(18a) 


where 


D i, j " ^ L i, 1 + R i, j + 


i,l 


+ B i,i } 


(18b) 


These L, R, T, B, and D expressions refer to the coefficients associated with 
any set of left, right, top, bottom and central mesh points as shown in figure 5. Be- 
cause of this point configuration, equation (18a) is called a five-point formula. This 


© 

u Uj 


© 

u i.i-l 


U U 


u® 
u i+l, j 


<D 

-- u i, j+i 

e 

Figure 5. - Five-point mesh configuration. 


equation is valid for all points except those on the polar axis. We will now discuss these 
exceptions by considering three special conditions; (1) the center point, (2) points on the 
9-0 positive polar axis, and (3) points on the 9 = n negative polar axis. 

(1) Center point (i = 1, 1 < j < m): In this case the cell simply becomes a sphere; 
and equation (18a) uses D- and R- coefficients only and becomes 

m m 

D i,i u i,i + Z R i,j u 2,rZ F i,i ' 1 (19a) 

i=i j=i 

(2) Points on the 0=0 axis (2 ^i s l- 1, j = 1): Setting 8q = 0j, equation (18a) 
reduces to 
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f 


(19b) 


D i, l u i, 1 + L i, l u i-l, 1 + R i, l u i+l, 1 + 2B i, l u i, 2 = F i, 1 


(3) Points on the 9 = tt axis (2 <i s 2-1, j = m): With = 0 m , equation (18a) 

reduces to 


D. „u, „ + L. __u- i __ + R. „u. , i „ + 2T- _u. _ - = F 


J' ___ u • _ -v j_i . u ■ h _ -r xv- u • , -t 

ij m i ? in. ij m i— 1 ? m ijin i+ 1 7 m 


i,m i,m- 


i,m 


(19c) 


Notice that in the last two cases, the cell becomes the frustum of a spherical cone. 

Summarizing then, for each interior mesh point (r. , 0.) where u. . is unknown, the 

1 J J 

finite difference expression is derived in terms of its four adjacent mesh points. Notice 
that three kinds of cells were used: 

(1) The curvilinear cube enclosing all mesh points interior to the sphere but not on 
the axis, 

(2) The sphere enclosing the center point, 

(3) The frustum of a spherical cone enclosing all axial points except the central 
point. 


Ordering of Mesh Points 

In this section we will explore how the ordering of mesh points affects the form of 
the coefficient matrix. Let us consider first case I, the full sphere configuration. 

There are two natural ways of point ordering, ring by ring or ray by ray. Both methods 
are shown below in figure 6 for a 35-point sample (5 r-rings and 7 0-rays). 



(a) Ring-by-ring ordering. 


(b) Ray-by-ray ordering. 


Figure 6. - Example of mesh point ordering. 
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The ring-by-ring ordering shown in figure 4 has considerable advantages. At mesh 
point 0, equation (19a) would use points 0, 1, 2, . . . . , 7. At mesh point 1, equation (19b) 
would use points 0, 8, 2, and 1 for the L, R, B, and D coefficients, respectively. Sim- 
ilarly, for mesh point 2, equation (18a) would use points 0, 9, 1, 3, and 2 for the L, R, 
T, B, and D coefficients, respectively. This method is then continued for points 3 to 6. 
For mesh point 7, equation (19c) would use points 0, 6, 7, and 14. This process is then 
continued clockwise and outward for each ring until the finite difference equation is 
written at all mesh points. The resultant coefficient matrix is shown in figure 7. Note 


0 

1 2 3 4 5 6 7 





0 

1 2 

8 




0 

1 2 3 

9 




0 

2 3 4 

10 




0 

3 4 5 

11 


0 


0 

4 5 6 

12 




0 

5 6 7 

13 




0 

6 7 

14 





1 

8 9 

15 




2 

8 9 10 

16 




3 

9 10 11 

17 




4 

10 11 12 

18 




5 

11 12 13 

19 




6 

12 13 14 

20 




7 

13 14 

21 





8 

15 16 

22 




9 

15 16 17 

23 




10 

16 17 18 

24 




11 

17 18 19 

25 




12 

18 19 20 

26 




13 

19 20 21 

27 




14 

20 21 

28 





15 

22 23 





16 

22 23 24 



0 


17 

23 24 25 





18 

24 25 26 

37 ) 




19 

25 26 27 

vf) 




20 

26 27 28 

(34) 




21 

27 28 

35 ) 


Figure 7. - Example of block tridiagonal coefficient matrix resulting from applying the 
five-point finite difference equation at each interior mesh point of the full sphere 
configuration. Note that the diagonal blocks are tridiagonal matrices, and the off- 
diagonal blocks are simply diagonal matrices. Circled numbers represent coeffi- 
cients of the surface boundary values. (The elements of any particular row are 
formed by L . . . TDB . . . R coefficients. For example, the elements of row 10 
(2 ... 8, 9, 10 ... 16) are formed by the left, top, central, bottom and right 
coefficients associated with mesh point 9 in fig. 4). 


that the partitioning of the rows and columns of the matrix in a 1, 7, 7, 7, 7 fashion 
results in nonzero elements lying in blocks that are zero except for the diagonal blocks, 
the superdiagonal blocks, and the subdiagonal blocks. We define this matrix as "block 
tridiagonal. " 

The other method of point ordering (ray-by-ray) shown in figure 6 would not yield 
a matrix with this block tridiagonal property. This is illustrated in figure 8. In a sim- 
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Figure 8. - Coefficient matrix resulting from improper choice of mesh point ordering. Notice the 
outriding coefficients in row 1 and column 1. Circled numbers represent coefficients of surface 
boundary values. 



Figure 9. - Example of ordering 
of mesh points for spherical 
segment (cases II and III) 
Boundary values are speci- 
fied on the spherical surface 
by points 29 to 35 and on the 
conical surface (0 Q = constant 
surface) by points 0, 1, 8, 15, 

22, and 29. Optional boundary 
values on the negative polar axis 
(0 = 7 r line) are specified by 
points 0, 7, 14, 21, 28, and 35. 



Figure 10. - Example of block tridiagonal coefficient matrix generated from the mesh 
points of the unknowns in the spherical segment configuration. Circled numbers repre- 
sent coefficients for the boundary values. 
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XDf 


2 3 


9 



2 3 4 


10 



3 4 5 


11 



4 5 6 


12 



5 6 

(JL 

13 



2 


9 10 


16 


3 


9 10 11 


17 


4 


10 11 12 


18 


5 


11 12 13 


19 


6 


12 13 

(g> 

20 



9 

© 

16 17 


23 



10 


16 17 18 


24 



11 


17 18 19 


23 



12 


18 19 20 


26 



13 


19 20 

<8L 

27 

(g 



16 

© 

23 24 

<56 



17 


23 24 25 

0 ] 



18 


24 25 26 

if? 



19 


25 26 27 

(33 



20 


26 27 



Figure 11. - Example of block tridiagonal coefficient matrix resulting from applying finite 
difference equation to mesh points of each unknown in example of spherical segment with 
axial values. Circled points represent coefficients for boundary values. 


ilar fashion, the point ordering and matrix generation for cases II and III are shown in 
figures 9 to 11. Note that the block tridiagonal property persists throughout, the only 
change being that of a decrease in the order of the matrices as one moves from case to 
case. 


Generation of Coefficient Matrix 

We have just shown with a small sample mesh how proper ordering of points is nec- 
essary for the formation of the block tridiagonal matrix. 

The corresponding coefficient matrix associated with case I, the full sphere config- 
uration, will now be examined. Here, ring-by-ring ordering of the 15 by 25 r-0 mesh 
points yields a block tridiagonal matrix whose nonzero elements all lie in a band of 
width 51 around the principal diagonal. This was shown previously in figure 7. 

It is possible at this time to write the matrix equation and, using an algorithm of 
the type given in reference 1, to obtain a solution. But this would necessitate working 
with a matrix requiring 19 125 storages (51 by 375). But, when a coefficient matrix 
exhibits the tridiagonal block property and, in addition, has off-diagonal blocks that are 
in themselves diagonal, then a modification of the algorithm enables one to solve the 
matrix equation at a considerable savings in storage. One needs to store only the in- 
verse of the block matrices associated with the principal diagonal blocks. For the pre- 
sent mesh size, this would require only 9375 storages (15 by 625). 
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Note also that the matrix is symmetrical. This is clearly shown by figure 7 and 
equation (15) where 


L i+lJ 



and 


T 


i, J+l 



From this property we may now write 


A 


i+1 



where the A's and C’s are block matrices as defined in equation (20). 

In this problem the five-point stencil previously shown in figure 5 produces a set of 
matrix equations MZ = K, where matrix M is block tridiagonal of the form 


and 



whose elements are matrices. Note that the A’s and C's are diagonal, the B's are 
tridiagonal symmetric, and K is simply a column vector. 


Development of Algorithm 

The matrix equation can be solved directly by using the algorithm of reference 1. 
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Matrices Wj, Wg, , W Q _i are calculated from 


W 1 - B I lc X 


S-l, 


(21a) 


W i = < B i - A i W i-l) C i 

Vectors Gj, Gg, , G n are calculated from 


2<i< (n-1) 


= b; 1 *! 


X-l 


(21b) 


G i " < B i ' A i W i-l> < K i - A i G i-l) 2 ^ i ^ 

Then the vector components Z. of the solution of MZ = K are given recursively by 

Z = G 
n n 


(21c) 


Z i = G i - W i Z i + l 


(n-1) > i > 1^ 

The algorithm may be modified as follows. Rewrite equation (21a) as 


<V A i w i-i>' 1 = w i c i' 1 


(22a) 


Now define 


v i = w i c i 


(22b) 


then 


W i-l= V i-l C i-l 


Since by symmetry A. = C-_ j, this may be rewritten as 


W i-1 = V i-l A i 


-1 


By definition = (B^ - A.Wj^) , then 


V i = < B i - A i V i-l A i) 


-1 


(22c) 


(22d) 


(22e) 
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Using equation (22e), the algorithm may now be rewritten as 


V 1 = 


Vi = (B t - A^.^) 


-1 


2 < i < n 


(23a) 


G 1 = Bj 1 K 1 (Unchanged) 
Ci^V-dq-A-G-.j) 


2 s i < n 




(23b) 


Z n = G r (Unchanged) 


Z i = G i - V i A i + l Z i + l 


(n-1) >121 


(23c) 


Notice that two multiplication operations have been eliminated, one in the calculation of 
Vp and one in the calculation of G^. For the full sphere, Bj is a scalar, and Ag is 
an NT by 1 column matrix (L- coefficients of innermost ring). The V- are full and 
symmetric, and these are all that must be stored (15 by 625 storages). 

The following section describes the FORTRAN program written from the algorithm 
just developed. 


THE FORTRAN PROGRAM 
General Description 

The program consists of nine subroutines: 

(1) POISS1 - generates r and 9 mesh arrays and sets all switches required by 

the spherical geometry specified. 

(2) SPHER1 - generates all coefficient matrices required for case I, the full 

sphere configuration. 

(3) SPHER2 - uses the algorithm to calculate the solution for case I. 

(4) SEG1 - generates all coefficient matrices required for either the case II 

or case m configuration 

(5) SEG2 - uses the algorithm to calculate the solution for either case II or 

case in. 

(6) FACTOR - performs matrix factorization before calculating inverses. 
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(7) INVERT - performs matrix inversions. 

(8) AE - calculates coefficients associated with each interior mesh point. 

(9) FE - calculates coefficients associated with each value of the source 

term array. 

The FORTRAN program is listed in appendix B and flow charts in appendix D. All 
FORTRAN symbols used are defined in appendix C. The program will return a solution 
in about 6 or 7 seconds and requires approximately 10 000 storages. 


Using the Program 

Using the subroutine is easy. First, the appropriate quantities are stored in 

AVS (logical constant) , NR, NT, XI, RA(NR), TA(NT), and X(NT, NR-1) 

Subroutine POISS1 is called. Then subroutine POISS2 is called, and the solution is re- 
turned in the X-array. This procedure will now be described in detail. 

Regardless of which case is being run, the main calling program must contain the 
following COMMON card 

COMMON/BLOCK1/NR, NT, RA(17), TA(25), X(25, 16), XI, AVS 
and the following definitions : 

NR number of r mesh lines, including the center point (maximum value, 17) 

NT number of 6 mesh lines (maximum value, 25) 

RA(IR) r mesh values, 1 s IR s NR (note that RA(1) = 0) 

TA(IT) 6 mesh values, 1 s IT s NT (note that TA(NT) = n) 

X(IT,NR-1) surface boundary values, 1 < IT :£ NT 

In addition, the main calling program must provide the following data: 

Case I, 9 - 0 (full sphere configuration): Load space charge terms into X(IT, IR), 

1 < IT < NT, 1 < IR < NR- 2 

Call POISS1 then POISS2. The solution is then returned in the X-array. 

Case II, 9 ± 0 (spherical cone; no axial values): 

(a) Load center point value into XI. 

(b) Load additional boundary values from the conical surface ( 9 = constant) into 
X(l, IR) (1 s IR < NR-2). 
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(c) Load space charge terms into X(IT, Ht) (2 s IT < NT; 1 < IR < NR-2). 

(d) Set AVS to false . 

Call POISS1, then POISS2. The solution is then returned in the X-array. 

Case in, 6*0 (spherical cone with axial values): 

(a) Load center point value into XI 

(b) Load additional boundary values from the conical surface (6 = constant) into 
X(l, IR) (1 < IR < NR-2). 

(c) Load axial boundary values into X(NT,IR) (1 s IR < NR-2). 

(d) Load space charge terms into X(IT,IR) (2 2 IT < NT-1, 1 <IR < NR-2). 

(e) Set AVS to true. 

Call POISS1, then POISS2. The solution is then returned in the X-array. Note that in 
all cases the solution is returned in the X-array, but the boundary values remain undis- 
turbed. When repeated solutions are sought for the same mesh but for different space 
charge terms, POISS1 need be called only once. 

A sample test run is given in appendix A which illustrates the running of each of the 
three cases described. 


CONCLUDING REMARKS 

We have described a numerical method and computer program for solving Poisson T s 
equation in axisymmetric regions of a sphere of two geometries: 

(1) A full sphere with boundary values on the surface only 

(2) A spherical cone with additional boundary values on the conical surface, and, 
optionally, along the negative polar axis. 

In developing the numerical method, a curvilinear cell was used at each mesh point 
to derive the finite difference analog to Poisson’s equation. The mesh points were then 
ordered in such a manner so as to produce a block tridiagonal symmetric coefficient 
matrix. The particular properties of this type of matrix were then used to write an 
algorithm from which a computer program was developed. The program uses a direct, 
noniterative solution and requires minimal storage of the requisite matrices. 

Two types of flexibility were built into the program: variable spherical geometry 
and irregular mesh spacing in two directions. 

It should be noted that the cell approach presents a simple way to avoid numerical 
difficulties in regions where singularities might be present. For example, in this prob- 
lem such a region would be along the polar axis. The numerical method is also general 
enough for easy extension to other equations and other geometries. 

The program is easy to use and general enough to be used as a black box in any 
situation where a solution to Poisson’s equation is needed in spherical regions. 
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For a mesh with a maximum allowable number of 376 points, the program returns a 
solution in 6 to 7 seconds (IBM DCS 7094/7044 computer) and requires an approximately 
10 000- word storage. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 14, 1971, 

129-04. 
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APPENDIX A 


SAMPLE PROBLEM 


A potential distribution was assumed to be of the form 


u = r cos 9 


(Al) 


Therefore, 




u = 3r cos 9 


u rr = 6r cos 9 


u^ = -r sin 9 


Uqq = -r cos 9 


(A2) 


Substituting these quantities into equation (2a), the expression for the right-hand 
side is obtained: 


6r cos 9 + — 3r 2 cos 9 + — (-r 3 cos 9) + - - ot — (-r 3 sin 0) = 10r cos 9 (A3) 


For all three cases, NR = 17, NT =25, and the RA array was 0, 0.2, 0.3, 0.39, 0.48, 
0.56, 0.63, 0.70, 0.76, 0.81, 0.86, 0.90, 0.93, 0.96, 0.98, 0.99, and 1.00. 

The TA array varied for all three cases in the following manner: 

Case I.(0 O = 0°) - TA = 0, 14, 28, 41, 54, 66, 77, 88, 98, 108, 117, 125, 133, 140, 
147, 153, 158, 163, 167, 171, 174, 176, 178, 179, 180. 

Case II (e Q = 45°) - TA = 45, 56, 66, 76, 86, 95, 103, 111, 119, 126, 133, 139, 145, 
150, 155, 160, 164, 168, 171, 174, 176, 177, 178, 179, 180. 

Case III (0 Q = 60°) - TA = 60, 70, 79, 88, 96, 104, 112, 119, 126, 132, 138, 144, 
149, 154, 158, 162, 166, 169, 172, 174, 176, 177, 178, 179, 180. 

Once these values are stored, the values on the right-hand side are calculated from 

equation (A3) and stored in X(IT, NR-1) (where IT = 1,2, ,25 and NR = 2, 3, .... , 

16). The following sets of boundary values are now calculated from equation (Al). For 
cases I, II, and m, the surface boundary values (r = 1) are calculated and stored in 
X(IT, 16) (where IT = 1,2,...., 25). For cases II and III, the conical surface 
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(i 6 = constant) boundary values are calculated and stored in X(l, IR) (where IR = 1, 2, . . 
16). For case IE only, the negative axial ( 9 = n) boundary values are calculated and 
stored in X(25, IR) (where IR = 1, 2, . . . . , 16). 

POISS1 is then called, followed by the final call, POISS2. The exact answers are 
calculated using equation (Al) and compared with those calculated by the subroutine. 
These residuals are also printed out for each case. Note that the largest residuals 
occur in those areas where the mesh spacing is coarsest. 

The FORTRAN program and output follow: 


SIBFTC TEST 

COMMON/ BLOCK 1 /NR , NT, R A l 1 /),TA(25),X(25,16),X1,AVS 
LOGICAL AVS 

NAME LIST/DATA/ TA.RA.NR.NT 
C 

14 FORMAT I 1X,F6.2, 16F7. 2 ) 

15 FORMAT ( 6H THETA , 3X, 1 6 I F5 .3 » 2X ) ) 

16 FORMAT! 1 1H1RES I DUALS- , 41 X, 3HRH0 ) 

17 FORMAT! 14H0AXIAL VALUES-) 

18 FORMAT 14X* 1HR , 8X , 1HV ) 

19 FORMAT l 2X»F5.3,3X»F6. 3) 

20 FORMAT (8H10UTPUT— »52X»3HKHQ) 

21 FORMAT! 14H0EXACT ANSWERS ,45X , 3HRH0 ) 

22 FORMAT! 1 X , F6 . 2 , 16F7. 4 ) 

23 FORMAT! 7H0 I NPUT-, 53X , 3HRH0) 

24 FORMAT! 1H1 ) 

25 FORMAT ( 55H ) 

26 FORMAT ! 6H XCTR=,F5.2) 

CVF=57. 29578 
DO 13 L = 1 » 3 
WRITE! 6,24) 

READ!5, 25 ) 

WRITE (6,25) 

REAl)(5»25 > 

WRITE! 6, 25) 

READ (5, DATA) 

AVS=. FALSE. 

NRMl^R-1 

CONVERT f HE T A ARRAY TO RADIANS. 

DO 1 J = 1 » NT 

1 T A ( J ) = T A ( Jl/CVF 

GENERATE RIGHT-HAND SIDE. 

DO 2 I T = 1 , NT 
TH = T A ( IT ) 

CT = C0S ( TH ) 

DO 2 IR=1,NRM1 

2 X(IT,IR)=10.*RA( IR+1 )*CT 
XI =0.3 
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C CALCULATE SURFACE BOUNDARY VALUES. 

DO 3 I T = 1 1 NT 

3 X( IT»16)=C0S!TA! IT)) 

IF(L.EQ.1 ) GO TO 7 

CALCULATE RADIAL BOUNDARY VALUES. 
CT=COS ( T A f 1 ) ) 

DO 4 I R = 1 » NRM 1 

4 X(l, IR ) =5 f RA f IR+1)**3)*CT 
IF(L.EQ.2) GO TO 7 

CALCULATE /PR I NT AXIAL BOUNDARY VALUES. 
AVS= • TRUE • 

CT = COS ! T A ( NT ) ) 

DO 5 I R= 1 » NRM 1 
R3= R A ( I R + 1 ) * * 3 

5 X ( NT * I R ) =R3*CT 
WRITE! 6,17) 

WR I T E ( 6 ? 1 8 ) 

WR I TF ! 6 , 1 9 ) R A ( 1 ) t XI 
DO 6 1=1* NRM 1 
I P 1 = I + 1 

6 WRITE! 6*19) R A ( I P 1 ) * X ( NT * I ) 

7 CALL PDISS1 

PRINT INPUT BLOCK. 

WRITE(6*23) 

WR I T E ( 6,15) ( R A ( I ) » 1=2, NR) 

DO 8 J = 1 , N T 
THD = CV F *T A { J ) 

8 WRITE(6*14) THD , (XtJ,K-l), K=2,NR) 

CALL POISS2 

PRINT OUTPUT BLOCK. 

WRI TE( 6,20 ) 

WR I TE ( 6 * 1 5 ) ( R A ( I ) , 1=2, NR) 

DO 9 J = 1 , NT 
THD = CVF*T A ( J ) 

9 WRITE(6,22) THD, ( X ( J , K- 1 ) * K = 2 , NR ) 

WR I T E ( 6 , 2 6 ) XI 

CALCULATE AND PRINT EXACT ANSWERS. 

WR I TE C 6 , 2 1 ) 

WRITE! 6,15) ! R A ( I ) , 1=2, NR) 

DO 10 1=1, NRM1 
DO 10 J = 1 , N T 
R 3= R A! I + 1 ) **3 
CT = COS t TA ( J ) ) 

RESIDUALS. 

DIMENSION RfcS ! 25, 16) 

T 1 = R3* C T 

RES! J, I )=X( J,I ) - T 1 
1C X ( J , I ) = T 1 

DO 11 J = 1 , N T 
THD=C V F *T A ( J ) 

11 WRITE(6,22) THD, !X(J,K-l), K=2,NR) 
WRITE!6,26) XI 
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C PRINT RESIDUALS. 

WRITE16.16) 

WRITE(6,I5) I R A l I ) » 1=2, NR) 

DO 12 J = 1 * NT 
THD=CVF*T A ( J ) 

12 WRI rE(6»22) THD, (RES{J,K-1), K=2,NR) 

13 CONTINUE 
C 

STOP 

END 
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TEST CASE l f THETA 0=0 DEGREES . 
(FULL SPHERE ). 


INPUT - 

THETA 

0.200 

0.300 

0.390 

0.480 

0.560 

0.630 

0.700 

RHO 

0.760 

0.810 

0.860 

0.900 

0.930 

0.960 

0.980 

0 . 99 J 

I . OjO 

0 . 

2.00 

3.00 

3.90 

4.80 

5.60 

6.30 

7.00 

7.60 

8.10 

8.60 

9.00 

9.30 

9.60 

9.80 

9.90 

1.00 

14.00 

1.94 

2.91 

3.78 

4.66 

5.43 

6.11 

6.79 

7.37 

7.86 

8.34 

8.73 

9.02 

9.31 

9.51 

9.61 

0.97 

28.00 

1.77 

2.65 

3.44 

4.24 

4.94 

5.56 

6.18 

6.71 

7.15 

7.59 

7.95 

8.21 

8.48 

8.65 

8 . /4 

0.88 

41.00 

1.51 

2.26 

2.94 

3.62 

4.23 

4.75 

5.28 

5.74 

6.11 

6.49 

6.79 

7.02 

7.25 

7.40 

7.47 

0.75 

54.00 

1.18 

1.76 

2.29 

2.82 

3.29 

3.70 

4.11 

4.47 

4.76 

5.05 

5.29 

5.47 

5.64 

5.76 

5. 82 

0.59 

66.00 

0.81 

1.22 

1.59 

1.95 

2.28 

2.56 

2.85 

3.09 

3.29 

3.50 

3.66 

3.78 

3.90 

3.99 

4.03 

0.41 

77.00 

0.45 

0.67 

0.88 

1.08 

1.26 

1.42 

L • 57 

1.71 

1.82 

1.93 

2.02 

2.09 

2.16 

2.20 

2.23 

0.22 

88.00 

0.07 

0.10 

0.14 

u .17 

0.20 

0.22 

0.24 

0.27 

0.28 

0.30 

0.31 

0.32 

0.34 

0.34 

0.35 

0.03 

98.00 

- 0.28 

- 0.42 

- 0.54 

- 0.67 

- 0.78 

- 0.88 

- 0.97 

- 1.06 

- 1.13 

- 1.20 

- 1.25 

- 1.29 

- 1.34 

- 1.36 

- 1 . 3 d 

- 0.14 

108.00 

- 0.62 

- 0.93 

- 1.21 

- 1.48 

- 1.73 

- 1.95 

- 2.16 

- 2.35 

- 2.53 

- 2.66 

- 2.78 

- 2.87 

- 2.97 

- 3.03 

- 3 . C 6 

- 0.31 

117.00 

- 0.91 

- 1.36 

- 1.77 

- 2.18 

- 2.54 

- 2. 86 

- 3.18 

- 3.45 

- 3.68 

- 3.90 

- 4.09 

- 4.22 

- 4.36 

- 4.45 

- 4.49 

- 0.45 

125.00 

- 1.15 

- 1.72 

- 2.24 

- 2.75 

- 3.21 

- 3.61 

- 4.02 

- 4.36 

- 4.65 

- 4.93 

- 5.16 

- 5.33 

- 5.51 

- 5.62 

- 5.68 

- 0.57 

133.00 

- 1.36 

- 2.05 

- 2.66 

- 3.27 

- 3.82 

- 4.30 

- 4.77 

- 5.18 

- 5.52 

- 5.87 

- 6.14 

- 6.34 

- 6.55 

- 6.68 

- 6.75 

- 0.68 

140.00 

- 1.53 

- 2.30 

- 2.99 

- 3.68 

- 4.29 

- 4.83 

- 5.36 

- 5.82 

- 6.23 

- 6.59 

- 6.89 

- 7.12 

- 7.35 

- 7.51 

- 7 . 58 

- 0 . 77 

147.00 

- 1.68 

- 2.52 

- 3.27 

- 4.03 

- 4.70 

- 5.28 

- 5.87 

- 6.37 

- 6.79 

- 7.21 

- 7.55 

- 7.80 

- 8.05 

- 8.22 

- 8.3 j 

- 0.84 

153.00 

- 1.78 

- 2.67 

- 3.47 

- 4.28 

- 4.99 

- 5.61 

- 6.24 

- 6.77 

- 7.22 

- 7.66 

- 8.02 

- 8.29 

- 8.55 

- 8.73 

- 8 . 82 

- U .89 

158.00 

- 1.85 

- 2.78 

- 3.62 

- 4.45 

- 5.19 

- 5.84 

- 6.49 

- 7.05 

- 7.51 

- 7.97 

- 8.34 

- 8.62 

- 8.90 

- 9.09 

- 9.18 

- 0.93 

163.00 

- 1.91 

- 2.87 

- 3.73 

- 4.59 

- 5.36 

- 6.02 

- 6.69 

- 7.27 

- 7.75 

- 8.22 

- 8.61 

- 8.89 

- 9.18 

- 9.37 

- 9 . 47 

- 0.96 

167.00 

- 1.95 

- 2.92 

- 3.80 

- 4.68 

- 5.46 

- 6.14 

- 6.82 

- 7.41 

- 7.89 

- 8.38 

- 8.77 

- 9.06 

- 9.35 

- 9.55 

- 9.65 

- C .97 

171.00 

- 1.98 

- 2.96 

- 3.85 

- 4.74 

- 5.53 

- 6.22 

- 6.91 

- 7.51 

- 8.00 

- 8.49 

- 8.89 

- 9.19 

- 9.48 

- 9.68 

- 9.78 

- 0.99 

174.00 

- 1.99 

- 2.98 

- 3.88 

- 4.77 

- 5.57 

- 6.27 

- 6.96 

- 7.56 

- 8.06 

- 8.55 

- 8.95 

- 9.25 

- 9.55 

- 9.75 

- 9. 85 

- 0.99 

176.00 

- 2.00 

- 2.99 

- 3.89 

- 4.79 

- 5.59 

- 6.28 

- 6.98 

- 7.58 

- 8.08 

- B .58 

- 8.98 

- 9.23 

- 9.58 

- 9.78 

- 9.88 

- 1.00 

178.00 

- 2.00 

- 3.00 

- 3.90 

- 4.80 

- 5.60 

- 6.30 

- 7.00 

- 7.60 

- 8.10 

- 8.59 

- 8.99 

- 9.29 

- 9.59 

- 9.79 

- 9.89 

- 1.00 

179.00 

- 2.00 

- 3.00 

- 3.90 

- 4.80 

- 5.60 

- 6.30 

- 7.00 

- 7.60 

- 8.10 

- 8.60 

- 9.00 

- 9.30 

- 9.60 

- 9.80 

- 9.90 

- 1 . 00 

180.00 

- 2.00 

- 3.00 

- 3.90 

- 4.80 

- 5.60 

- 6.30 

- 7.00 

- 7.60 

- 8.10 

- 8.60 

- 9.00 

- 9.30 

- 9.60 

- 9.80 

- 9.90 

- l.uO 


to 

CD 


03 

O OUTPUT- RHO 

THETA 0.203 0.300 0.390 0.480 0.560 0.630 0.700 0.760 0.810 0.860 0.900 0.930 0.960 0.980 C.990 1.000 

0. 0.0089 0.0282 0.0605 0.1114 0.1762 0.2504 0.3432 0.4390 0.5314 0.6360 0.7289 0.8043 0.8847 0.9412 0.9703 1.0000 

14.00 0.0086 0.0274 0.0587 0. 1081 0.1709 0.2430 0.3329 0.4260 0.5156 0.6171 0.7073 0.7804 0.8584 0.9132 0.9415 0.9703 

28.00 0. 007 B 0.0249 0.0534 0.0983 0.1555 0.2210 0.3029 0.3875 0.4692 0.5615 0.6435 0.7101 0.7811 0.8310 0.8567 0.8829 

41.00 0.0066 0.0212 0.0456 0.0840 0.1328 0.1889 0.2589 0.3312 0.4010 0.4799 0.5501 0.6070 0.6677 0.7103 0.7323 0.7547 

54.00 0.0051 0.0164 0.0354 0.0653 0.1033 0.1470 0.2015 0.2578 0.3122 0.3737 0.4283 0.4727 0.5200 0.5532 0.5703 0.5878 

66.00 0.0034 0.0113 0.0244 0.0451 0.0714 0.1016 0.1393 0.1783 0.2160 0.2585 0.2963 0.3270 0.3598 0.3828 0.3946 0.4067 

77.00 0.0018 0.0061 0.0133 0.U248 0.0394 0.0561 0.077C 0.0986 0.1194 0.1429 0.1639 0.1809 0.1990 0.2117 0.2183 0.2250 

88.00 0.0000 0.0007 0.0018 0.0036 0.0058 0.C085 0.0117 0.0151 0.0183 0.0220 0.0253 0.0280 0.0308 0.0328 0.0338 0.0349 

98.00-0.0016-0. 0043-0.008 8-0.0158-0.0248-0.0351-0.0480-0.0613-0 .0741-0.0886-0. 1015-0. 1120-0. 1232 -0.131 0-0. 13 50—0 .1392 

108.00- 0. 003 1-0 .0091 -0.0191 -0.0348-0. 0548-0. 0777-0. 106 3- 0.1359-0. 1644-0.1 967-0.2254-0. 2486-0.2 734-0. 2909-0 .2998-0* 3090 

117. 00- 0. 0045-0. 0132-0. 0278-0. C509-0. 0803-0. 1140-0. 1560-0. 1995-0. 2414-0. 2 889-0. 3310-0. 3652-0. 4 017-0. 4273- 0.44C5-0. 4540 

125.00- 0.0056-0.0166-0.0351-0.0642-0.1013-0.1439-0.1970-0.2520-0.3049-0.3649-0.4181-0.4614-0.5075-0.5398-0.5565-0.5736 

133.00- 0.0065-0.0197-0.0416-0.0763-0.1204-0.1710-0.2342-0.2995-0.3625-0.4338-0.4972-0.5486-0.6034-0.6419-0.6617-0.6820 

140.00- 0.0073-0.0221-0.0467-0.0856-0.1352-0.1920-0.2630-0.3364-0.4072-0.4872-0.5584-0.6161-0.6777-0.7210-0.7433-0.7660 

147.00- 0.0080-0.0241-0.0511-0.0937-0.1480-0.2102-0.2879-0.3683-0.4458-0.5334-0.6113-0.6746-0.7420-0.7893-0.8130-0.8387 

153.00- 0.0084-0.0256-0.0543-0.0995— 0.1572^-0. 22 33-0. 3058-0. 39 12-0 .473 6- 0.5667- 0.649 5-0.716 6-0.7 88 3-0. 8386-0 .8645—0 .8910 

158.00- 0.0088-0.0266-0.0565-0. 1036-0. 1635-0. 2323-0.3 182-0. 4071-0. 4928-0. 5897-0. 6758-0. 7457-0. 8203-0. 8726-0 .89 S6-0 .9272 

163.00- 0.0090-0.0274-0.0582-0. 1068-0.1686-0.2396-0.3282-0.4198-0.5082-0.6082-0.6971-0.7691-0.8460-0.9000-0.9279-0.9563 

167.00- 0. 0092-0.0279-0.0593-0. 1088-0.1718-0.2441-0.3344-0.4278-0.5178-0.6197-0.7102-0.7837-0.8620-0.9170-0.9454-0.9744 

171. 00- 0. 0093-0. 0283-C. 060 1-0. 1103-0.1742-0.2475-0.3390-0.4336-0.5249-0.6281-0.7199-0.7944-0.8738-0.9296-0.9583-0.9877 

174.00- 0.0094-0.0285-0.0605-0. 1110-0.1754-0.2492-0.3413-0.4366-0.5285-0.6325-0.7249-0.7999-0.8798-0.9360-0.9650-0.9945 

176.00- 0.0094-0.02 86-0.0607-0. 11 14-0. 1759-0. 2499-0.3423-0.4379-0. 5 301-0. 6 344-0. 727 1-0. 80 2 3-0. 8825-0. 9389-0 .96 79-0 . 99 76 

178.00- 0.0094-0.0286-0.0608-0. 1116-0.1762-0.2504-0.3430-0.4387-0. 531 1- 0 . 635 6-0. 7284-0 . 0038-0. 8 841 -0. 9406-C .9697-0 . 9994 

179.00- 0.0094-0.0287-0.0608-0. 1116-0.1763-0.2505-0.3431-0.4389-0.5313-0.6358-0. 7288-0. 8041-0. 8 845-0. 9410-0. 970 1-0 . 9998 

180.00- 0.0094-0.0287-0.0608-0. 11 16-0 . 1 763-0. 2505-0 . 3432-0. 4390-0 . 53 14-0 . 6359-0 . 7289-0. 8D43-0. 8 847-0. 941 2-0 .97 G3-1 . 0000 
XCTR=-0. 00 


EXACT ANSWERS RHO 


THETA 

0.203 

0.300 

0.390 

0.480 

0.560 

0.630 

0.700 

0.760 

0.810 

0.860 

0.900 

0.930 

0.960 

0.980 

C.990 

l.OJO 

0 . 

0.0083 

0.0270 

0.0593 

0.1106 

0.1756 

0.2500 

0.3430 

0.4390 

0.5314 

0.6361 

0.7290 

0.8044 

0.8847 

0.9412 

0.9/C3 

1.0000 

14.00 

0.0078 

0.0262 

0.0576 

0. 1073 

0.1704 

0.2426 

0.3328 

0.4259 

0.5157 

0.6172 

0.7073 

0.7805 

0.8585 

0.9132 

0.9415 

0.9703 

28.00 

0.0071 

0.0238 

0.0524 

0.0976 

0.1551 

0.2208 

0.3029 

0.3876 

0.4692 

0.5616 

0.6437 

0.7102 

0.7812 

0.8310 

C • 85 67 

0 . 88 29 

41.00 

0.0060 

0.0204 

0.0448 

0.0835 

0.1325 

0.1887 

0.2589 

0.3313 

0.4011 

0.4800 

0.5502 

0.6071 

0.6677 

0.7103 

0.7^23 

0.7547 

54.00 

0.0047 

0.0159 

0.0349 

0.0650 

0.1032 

0.1470 

0.2016 

0.2580 

0.3124 

0.3739 

0.4285 

0.4728 

0.5200 

0.5532 

0.5703 

0.58 78 

66.00 

0.0033 

0.0110 

0.0241 

0.C450 

0.3714 

0.1017 

0.1395 

0.1785 

0.2162 

0.2587 

0.2965 

*.3272 

0.3 599 

0.3828 

0.3947 

0.4067 

77.00 

0.0018 

0.0061 

0.0133 

0 • 0249 

0.U395 

0.0562 

0.0772 

0.0987 

0.1195 

0.1431 

0.1640 

0.1809 

0.1990 

0.2117 

0.2183 

0.2250 

88.00 

0.0003 

0.0009 

0.0021 

0.0039 

0.0061 

0.0387 

0.0120 

0.0153 

0.0185 

0.3222 

0.0254 

0.0281 

0.0309 

0.0328 

0.0339 

0.0349 

98.00- 

-0.0011- 

-0 .00 38- 

-0.0083- 

-0.0154- 

-0.0244- 

-0.0348- 

-0.0477- 

-0.0611- 

-0.0740- 

-0.0885- 

-0.1015- 

-0.1119- 

-0.1231- 

-0.1310- 

-0.135C- 

-0.1392 


108.00- 0*0025-0.0083-0*0183-0 342-0. 0543-0. 0773-0. 1060-0. 1357-0. 1642-0. 196 6- 0.2253-0. 2^86-0. 2734-0.2908-0.2998-0. 3090 

117.00- 0.0036-0. Cl 23-0. 0269-0. 0502-0.0797-0. 1135-0.1557-0.1993-0.2413-0.2888-0.3310-0.3652-0.4017-0.4273-0.4405-0.4540 

125.00- 0.00^6-0.0155-0.0340-0. 0634-0.1007-0. 1434-0.1967-0.2518-0.3048-0.3648-0.4181-0.4614-0.5075-0.5398-0.5565-0.5736 

133.00- 0.0055-0.01 84-0.0405-0. 0754-0.1198-0. 1705-0. 2339-0. 29 94-0. 362 4- 0.^338-0 .497 2-0. 5486-0. 6034-0 .6 41 9- 0.661/- 0.68 20 

140 .00- 0 .006 1-0 .02 07-0 .045 4-0. 0847-0. 1345-0. 1915-0.2628-0.3363-0. 4071-0.^872-0. 5584-0. 6162-0. 6777-0. 7210-C .7433-0 .7660 

147.00- 0.0067-0.0226-0.0497-0.0928-0. 1473-0.2097-0.2877-0.3682-0.4457-0.5334-0.6114-0.6746-0.7420-0.7893-0.8138-0.8387 

153. 00- 0. 307 1-0. 0241-0. 052 9-0. 09 8 5-0. 1565-0. 2228-0. 3056-0. 391 1-0. 4735-0. 5 667-0. 6495-0. 71 67-0. 7 883-0. 8386-0. 8645-0. 89 10 

158.00- 0. 007^-0.0250-0.0550-0. 1025-0 . 1628-0 . 23 18-0 . 3 180-0 . 4070-G . 4927-0 . 5 897-0 . 6759-0 . 7458-0 . 8 203-0. 8727- 0 . 89 96-0 . 92 72 

163.00- 0.0077-0.0258-0.0567-0. 1058-0.1679-0.2391-0.3280-0.4198-0.5082-0.6083-0.6971-0.7692-0.8461-0.9001-0.9279-0.9563 

167.00 -0. 0 07 8-0. 0263-0. 057 8-0. 1078-0.1711-0.2436-0.3342-0.4277-0. 51 78-0. 6 198-0. 7 103-0 . 7837-0. 8 62 1 -0 . 917 1-0 . 9454-0 . 97 44 

171.00- 0. D079-0 .0267- 0.0586 -0.1092-0. 1735-0. 2470-0. 3388-0. 43 3 6-0. 5249- 0.6282-0.720 0-0.7945-0. 8 73 8 -0.929 6-0.95 84-0 .9877 

1 74. 00- 0. 0C8D-0. 02 69-0. 0590-0. 1100-0.1747-0.2487-0.3411-0.4366-0.5285-3.6326-0.7250-0.8000-0.8799-0.9360-0 .9650-0.9945 

176.00- 0.0080-0.0269-0.0592-0. 1103-0.1752-0.2494-0.3422-0.4379-0.5301-0.6345-0.7272-0.8024-0.8826-0.9389-0.9679-0.9976 

178.00- 0.0083-0.0270-0.0593-0. 1105-0.1755-0.2499-0.3428-0.4387-0.5311-0.6357-0.7286-0.8039-0.8842-0.9406-0.9697-0.9994 

179.00- 0.0080-0.0270-0.0593-0. 1106-0 . 1756-0. 2500-0 . 3429-0 . 4389-0 . 5314- 0. 6360-0. 7289-0. 8042-0. 8 846-0 . 9410-0 .97C2-0 . 9998 

1 80.00- 0. 0080-0 .0270-0.0593-0. 1106-0 • 1756-0. 2 5C0-0. 3430-0. 4390-0. 5314-3. 6361-0. 7290-0. 8344-0. 8 847-0.9412-0 .9 / ^ 3-1 .0000 
XCTR=— 0 .00 


RESIDUALS- RHO 

THETA 0.203 0.300 0.390 0.480 0.560 0.630 0.700 0.760 0.810 0.860 0.900 0.930 0.960 0.980 0.990 1.000 

0. 0.0009 0.0012 0.0012 0.0008 0.0005 0.0004 0.0002 0.0000 0.0000-0.0001-0.0001-0.0001-0.0000-0.0000-0.0000 0. 

14.00 0.0009 0.0012 0.0011 0.0008 0.0005 0.0004 0.0001 0.0000-0.0000-0.0001-0.0001-0.0001-0.0000-0.0000-0.0000 0. 

28.00 0.0008 0.0010 0.0010 0.0006 0.0004 0.0003 0.0000-0.0001-0.0001-0.0001-0. 0001-0. 0001-0.0001-0. 0000-0 .0000 0. 

41.00 0.3006 0.0008 0.0008 0.0005 0.0003 0.0002-0. 0000-0. 0001-0. 0001-0.0001-0. 0001-0.0001-0.0001-0. 0000-0. OOC.O 0. 

54.00 0.0004 0.0006 0.0005 0.0003 0.0001 0.0000-0. 0001-0. 0002-0. 0002-0. 0002-0. 0002-0. 0001-0. 0001-0. 0000-0. OOCO 0. 

66.00 0.0002 0.0003 0.0002 0. 0001-0.0000-0.0001-0 .0002-0.0002-0.0002-0. 0002-0.0002-0. 0001-0.0001-0. 0000-0 .OOCO 0. 

77.00- 0.0000 0.0000-0. 0000-0. 0001-0. 0001-0. 0002-0. 0002-0. 0002-0. 0002-0. 0001-0. 000 1-0. 0001-0. 0000-0. 0000-0. OCcO 0. 

88. 00- 0. 0003-0. 0003-0. 0003-0. 0003-0. 0003-0. 0003-0. 0003-0. 0002-0. 0002-0. 0002-0. 0001-0. 0001-0. 0001 -0.0000-0. 00 CO 0. 

98.00- 0.0005-0.0005-0.0005-0.0004-0.0004-0.0003-0.0003-0.0002-0.0002-0.3001-0.0001-0.0001-0.0000-0.0000-0.0000-0. 

108.00- 0.0007-0.0008-0.0007-0.0006-0.0005-0.0004-0.0003-0.0002-0.0002-0.0001-0.0001-0.0001-0.0000-0.0000-0.0 CO-O. 

117. 00- 0. 0008-0. 00 10-0. 0009-0. 0007-0. 0006-0. 0005-0. 0003-0. 0002-0. 0002-0. 0001-0. 0001-0. 0000-0. 0000 -0.0000-0. 0000-0. 

125. 00- 0. 001 0-0. 001 1-0. 001 1-0. C008-0. 0006-0. 0005-0. 0003-0. 0002-0. 0001-0. 0000-0. 0000-0. 0000 0.0000 0.0000 O.OOCJ-O. 

133.00- 0.0011-0.0013-0.0012-0.0009-0.0006-0.0005-0.0003-0.0002-0.0001-0.0000-0.0000 0.0000 0.0000 0.0000 C.OOCO-O. 

140.00- 0.0012-0.0014-0.0013-0.0009-0.0007-0.0005-0.0003-0.0001-0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 O.OOCO-O. 

147.00- 0.0013-0.0015-0.0014-0.0010-0.0007-0.0005-0.0002-0.0001-0.0001 0.0000 0.0003 0.0000 0.0000 0.0000 0.0000-0. 

153.00- 0.0013-0.0015-0.0014-0.0010-0.0007-0.0005-0.0002-0.0001-0.0000 0.0000 0.0001 0.0000 0.0000 0.0000 0.0000-0. 

158.00- 0.0014-0.0016-0.0015-0.0010-0.0007-0.0005-0.0002-0.0001-0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0.0 CO-O. 

163. 00- 0. 0014-0. 0016-0. 0015-0. 0010-0. 0007-0. 0005-0. 0002-0. 0001 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0.0000-0. 

167.00- 0.0014-0.0016-0.0015-0.0010-0.0007-0.0005-0.0002-0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 O.OOCO-O. 

171.00- 0.0014-0.0017-0.0015-0.0010-0.0007-0.0005-0.0002-0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 O.O.CO-O. 

174. 00- 0. 0014-0. 0017-0. 0015-0. C010-0. 0007-0. 0005-0. 0002-0. 0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 O.OOCO-O. 

176.00- 0.0014-0.0017-0.0015-0. G010-0. 0007-0. 0005-0. 0002-0. 0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 O.OOCC-O. 

178. 00- 0. 0014-0. 0017-0. 0015-0. C011-0. 0007-0. 0005-0. 0002-0. 0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 O.OOCO-O. 

179.00- 0.0014-0.0017-0.0015-0.0011-0.0007-0.0005-0.0002-0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0.0000-0. 

180.00- 0.0014-0.0017-0.0015-0.0011-0.0007-0.0005-0.0002-0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 O.OOCC-O. 


co 



oo 

to 


TEST CASE 2 , THET AQ =45 DEGREES . 
(NO AXIAL BOUNDARY VALUES ). 


INPUT - 

THETA 

0.200 

0.300 

0.390 

0.480 

0.560 

0.630 

0.700 

RHO 

0.760 

0.810 

0.860 

0.900 

0.930 

0.960 

0.980 

0.990 

1.000 

45.00 

0.01 

0.02 

0.04 

0.08 

0.12 

0.18 

0.24 

0.31 

0.38 

0.45 

0.52 

0.57 

0.63 

0.67 

0.69 

0 . 71 

56.00 

1.12 

1.68 

2.18 

2.68 

3.13 

3.52 

3.91 

4.25 

4.53 

<*.81 

5.03 

5.20 

5.37 

5.48 

5.54 

0 . 56 

66.00 

0.81 

1.22 

1.59 

1.95 

2.28 

2.56 

2.85 

3.09 

3.29 

3.50 

3.66 

3.78 

3.90 

3.99 

4.03 

0.41 

76.00 

0.48 

0.73 

0.94 

1.16 

1.35 

1.52 

1.69 

1.84 

1.96 

2.08 

2.18 

2.25 

2.32 

2.37 

2.40 

0.24 

86.00 

0.14 

0.21 

0.27 

0.33 

0.39 

0.44 

0.49 

0.53 

0.57 

0.60 

0.63 

0.65 

0.67 

0.68 

0.69 

0.07 

95.00 

- 0.17 

- 0.26 

- 0.34 

- 0.42 

- 0.49 

- 0.55 

- 0.61 

- 0.66 

- 0.71 

- 0.75 

- 0.78 

- 0.81 

-0 .84 

- 0.85 

- 0.86 

- 0.09 

103.00 

- 0.45 

- 0.67 

- 0.88 

- 1.08 

- 1.26 

- 1.42 

- 1.57 

- 1.71 

- 1.82 

- 1.93 

- 2.02 

- 2.09 

- 2.16 

- 2.20 

- 2.23 

- 0.22 

HI . 00 

- 0.72 

- 1.08 

- 1.40 

- 1.72 

- 2.01 

- 2.26 

- 2.51 

- 2.72 

- 2.90 

- 3.08 

- 3.23 

- 3.33 

- 3.44 

- 3.51 

- 3.55 

- 0.36 

119.00 

- 0.97 

- 1.45 

- 1.89 

- 2.33 

- 2.71 

- 3.05 

- 3.39 

- 3.68 

- 3.93 

- 4.17 

- 4.36 

- 4.51 

- 4.65 

- 4.75 

- 4.80 

- 0 . 48 

126.00 

- 1.18 

- 1.76 

- 2.29 

- 2.82 

- 3.29 

- 3.70 

- 4.11 

- 4.47 

- 4.76 

- 5.05 

- 5.29 

- 5.47 

- 5.64 

- 5.76 

- 5.82 

- 0.59 

133.00 

- 1.36 

- 2.05 

- 2.66 

- 3.27 

- 3.82 

- 4.30 

- 4.77 

- 5.18 

- 5.52 

- 5.87 

- 6.14 

- 6.34 

- 6.55 

- 6.68 

- 6.75 

- 0.68 

139.00 

- 1.51 

- 2.26 

- 2.94 

- 3.62 

- 4.23 

- 4.75 

- 5.28 

- 5.74 

- 6 .il 

- 6.49 

- 6.79 

- 7.02 

- 7.25 

- 7.40 

- 7.47 

- 0.75 

145.00 

- 1.64 

- 2.46 

- 3.19 

- 3.93 

- 4.59 

- 5.16 

- 5.73 

- 6.23 

- 6.64 

- 7.04 

- 7.37 

- 7.62 

- 7.86 

- 8.03 

- 8.11 

- 0.82 

150.00 

- 1.73 

- 2.60 

- 3.38 

- 4.16 

- 4.85 

- 5.46 

- 6.06 

- 6.58 

- 7.01 

- 7.45 

- 7.79 

- 8.05 

- 8.31 

- 8.49 

- 8 . 57 

- 0.87 

155.00 

- 1.81 

- 2.72 

- 3.53 

- 4.35 

- 5.08 

- 5.71 

- 6.34 

- 6.89 

- 7.34 

- 7.79 

- 8.16 

- 8.43 

- 8.70 

- 8.88 

- 8.97 

- 0.91 

160.00 

- 1.88 

- 2.82 

- 3.66 

- 4.51 

- 5.26 

- 5.92 

- 6.58 

- 7.14 

- 7.61 

- 8.08 

- 8.46 

- 8.74 

- 9.02 

- 9.21 

- 9.30 

- 0.94 

164.00 

- 1.92 

- 2.88 

- 3.75 

- 4.61 

- 5.38 

- 6.06 

- 6.73 

- 7.31 

- 7.79 

- 8.27 

- 8.65 

- 8.94 

- 9.23 

- 9.42 

- 9.52 

- 0.96 

168.00 

- 1.96 

- 2.93 

- 3.81 

- 4.70 

- 5.48 

- 6.16 

- 6.85 

- 7.43 

- 7.92 

- 8.41 

- 8.80 

- 9.10 

- 9.39 

- 9.59 

- 9.68 

- 0.98 

171.00 

- 1.98 

- 2.96 

- 3.85 

- 4.74 

- 5.53 

- 6.22 

- 6.91 

- 7.51 

- 8.00 

- 8.49 

- 8.89 

- 9.19 

- 9.48 

- 9.68 

- 9.78 

- 0 . 99 

174.00 

- 1.99 

- 2.98 

- 3.88 

- 4.77 

- 5.57 

- 6.27 

- 6.96 

- 7.56 

- 8.06 

- 8.55 

- 8.95 

- 9.25 

- 9.55 

- 9.75 

- 9.85 

- 0.99 

176.00 

- 2.00 

- 2.99 

- 3.89 

- 4.79 

- 5.59 

- 6.28 

- 6.98 

- 7.58 

- 8.08 

- 8.58 

- 8.98 

- 9.28 

- 9.58 

- 9.78 

- 9.88 

- 1.00 

177.00 

- 2.00 

- 3.00 

- 3.89 

- 4.79 

- 5.59 

- 6.29 

- 6.99 

- 7.59 

- 8.09 

- 8.59 

- 8.99 

- 9.29 

- 9.59 

- 9.79 

- 9.89 

- rl.OO 

178.00 

- 2.00 

- 3.00 

- 3.90 

- 4.80 

- 5.60 

- 6.30 

- 7.00 

- 7.60 

- 8.10 

- 8.59 

- 8.99 

- 9.29 

- 9.59 

- 9.79 

- 9.8 9 

- 1.00 

179.00 

- 2.00 

- 3.00 

- 3.90 

- 4.80 

- 5.60 

- 6.30 

- 7.00 

- 7.60 

- 8.10 

- 8.60 

- 9.00 

- 9.30 

- 9.60 

- 9.80 

- 9.90 

- 1.00 

180.00 

- 2.00 

- 3.00 

- 3.90 

- 4.80 

- 5.60 

- 6.30 

- 7.00 

- 7.60 

- 8.10 

- 8.60 

- 9.00 

- 9.30 

- 9.60 

- 9.80 

-9 . 9 j 

- 1.00 



OUTPUT- RHO 

THETA 0.200 0.300 0.390 0.480 0.560 0.630 0. 700 0.760 0.810 0.860 0.900 0.930 0.960 0.980 C. 990 1.000 

45.00 0.0057 0.0191 0.0419 0.0782 0.1242 0.1768 0.2425 0.3104 0.3758 0.4498 0.5155 0.5688 0.6256 0.6655 0.6861 0.7071 

56.00 0.0044 0.0151 0.0332 0.0618 0.0981 0.1397 0.1917 0.2453 0.2970 C.3555 0.4075 0.4497 0.4947 0.5263 0.5 26 0.5592 

66.00 0.0031 0.0109 0.0241 0.0449 0.0713 0.1016 0.1393 0.1784 0.2160 0.2586 0.2964 0.3271 D.3598 0.3828 C.3946 0.4067 

76.00 0.0016 0.0063 0.0142 0.0265 0.0423 0.0603 0.0828 0.1060 0.1284 0.1537 0.1763 0.1945 0.2140 0.2277 0.2347 0.2419 

86.00 0.0001 0.0015 0.0038 0.0074 0.0119 0.0171 0.0236 0.0304 0.0369 0.0442 0.0507 0.0560 0.0617 0.0656 0.0677 0.0698 

95.00- 0.001 3-0.0029-0.0057—0. 0101 -0.0157-0. 0221 -0.0 30 2-0.038 5- 0.0 46 5- 0.055 6- 0.063 7- 0.0702 -0.0 772 -0.0 821-0 .0846-0.0872 

103.00- 0.0025-0.0068-0.0140-0.0254-0.0399-0.0566-0.0774-0.0989-0.1197-0.1432-0. 1640-0. 1810-0. 1990-0.2117-0.2183-0.2250 

111. 00- 0. 0037-0. 0106-0. 022 1-0. 0403-0. 0634-0. 0900-0. 1232-0. 15 75-0. 190 6- 0.2 28 0-0. 261 3-0. 28 8 3-0. 3 171 -0.337 3-0. 347 7-0. 35 84 

119. 00- 0. 0048-0. 0141-0. 0297-0. 0543-0. 0857-0. 12 16-0. 1665-0. 2 130-0. 2578- 0.3 084-0. 3 53 5-0. 3900-0. 4 289-0. 4563-0. 47 -:4-0. 48 48 

126.00- 0.0058-0.0170-0.0360-0.0658-0.1038-0.1474-0.2019-0.2582-0.3125-0.3739-0.4285-0.4728-0.5200-0.5532-0.57:3-0.5878 

133.00- 0.0066-0.0197-0.041 6-0. C763-0. 1204-0. 1710-0 .2342-0 . 2995-0. 3625-0. 4338-0. 4972-0. 5486-0. 6034-0. 6419-0.661 7-0.6820 

139. 00- 0. 0073-0. 0218-0. 0460-G. C844-0. 1332-0. 1892-0.2591-0.3314-0.401 1-0.4800-0. 5501-0.6070-0.6677-0.7103-0.7323-0.7547 

145.00- 0.007 3-0.0236-0.0499-0.0915-0.1445-0.2053-0.2812-0.3597-0.4354-0.5 210-0.5971-0.658 8-0.7 247-0.7710-0 .79 48-0.8192 

150. 00- 0. 0082-0. 0249-0. 0528-0. 0967-0. 1527-0. 2 170-0. 2972-0. 3802-0. 4602-0. 5 508-0. 63 12-0. 69 65-0. 7 662-0. 815 1-0. 8 (03-0.86 60 

155.00- 0.0086-0.0260-0.0552-0. 1012-0. 1598-0.2271-0. 3110-0.3979-0. 481 6-0.5 764-0. 6606-0. 72 89-0.8 018-0. 8530-0 .8 794-0. 9063 

160.00- 0.0009-0.0270-0.0572-0. 1049-0 . 165 7-0. 2354-0 . 3225-0 . 4125-0. 4994- 0. 5 976-0. 6849-0. 75 58-0. 8313-0. 8844-C .91 1 8-0 . 93 97 

164.00- 0.0091-0.0276-0.0585-0. 1073-0 . 1695-0. 2408-0 . 3299-0. 4220-0 . 5108-0. 6 1 1 3-0. 7006-0. 77 31-0 . 8 504-0. 9047-0 .93 27-0 . 96 13 

168.00- 0.0093-0.0281-0.0595-0. 1092-0 . 1724-0. 2450-0 . 3357-0. 4294-0. 5 198-0 . 6220-0. 71 29-0. 7867-0. 8653-0. 9206-0 .9491-0 . 9781 

171.00- 0.0093-0.0283-0.0601-0.1103-0.1741-0.2 474-0.3389-0.4336-0. 52 49-0.6 281-0. 7199-0.7944-0.8 738-0.9296-0.958 3-0.9877 

174.00- 0.0094-0.0285-0.0605-0. 11 10-0. 1753-0.2491-0. 3413-0.4366-0. 5285-0.6324-0. 7249-0.7999-0.8798-0.9360-0.9650-0.9945 

176.00- 0.0094-0.0286-0.0607-0. 11 13-0. 1759-0. 2499-0 . 3423-0 . 4379-0 . 5 30 1-0. 6344-0. 727 1-0. 8023-0. 8 825-0. 9389-0 .96 79-0 . 9976 

177.00- 0.0094-0.0286-0.0608-0. 1115-0. 1760-0.2502-0. 3427-0.4384-0. 5307-0. 635 1-0. 7279-0. 8032-0. 8834-0. 9399-C .9689-0.9986 

178.00- 0.0094-0.0287-0.0608-0. 1116-0.1762-0.2503-0.3429-0.4387-0.5311-0.6355-0. 7284-0.8038-0.8841-0.9406-0.9697-0.9994 

179.00- 0.0094-0.0287-0.0608-0. 1116-0.1763-0.2505-0.3431-0.4389-0.5313-0.6358-0.7288-0.8041-0.8845-0.9410-0.9701-0.9998 

180.00- 0.0095-0.0287-0.0608-0. 1116-0.1763-0.2505-0.3432-0.4390-0.5314-0.6359-0.7289-0.8043-0.8847-0.9411-0 .9703-1.0000 
XCTR= 0. 

EXACT ANSWERS RHO 

THETA 0.200 0.300 0.390 0.480 0.560 0.630 0.700 0.760 0.810 0.860 0.900 0.930 0.960 0.980 0.990 1.000 

45.00 0.0057 0.0191 0.0419 0.0782 0.1242 0.1768 0.2425 0.3104 0.3758 0.4498 0.5155 0.5688 0.6256 0.6655 0.6H61 0.7071 

56.00 0.0045 0.0151 0.0332 0.0618 0.0982 0.1398 0.1918 0.2455 0.2972 0.3557 0.4077 0.4498 0.4947 0.5263 0.5426 0.5592 

66.00 0.0033 0.0110 0.0241 0.0450 0.0714 0.1017 0.1395 0.1785 0.2162 0.2587 0.2965 0.3272 0.3599 0.3828 0.3947 0.4067 

76.00 0.0019 0.0065 0.0144 0.0268 0.0425 0.0605 0.0830 0.1062 0.1286 0.1539 0.1764 0.1946 0.2140 0.2277 0.2347 0.2419 

86.00 0.0006 0.0019 0.0041 0.0077 0.0123 0.0174 0.0239 0.0306 0.0371 0.D444 0.0509 0.0561 0.0617 0.0657 0.0677 0.0698 

95. 00- 0. O0D7-O. 0024-0. 0052-0. 0096-0. 0153-0. 02 1 8-0. 0299-0. 0383-0.0 463-0. D554-0. 0635-0. 0701-0. 0 771-0. C820-0. 084 6-0. 08 72 

103.00- 0.0018-0.0061-0.0133-0.0249-0.0395-0.0562-0.0772-0.0987-0.1195-0.1431-0.1640-0.1809-0.1990-0.2117-0.2183-0.2250 

111. 00- 0. 0029-0. 0097-0. 021 3-0. 0396-0. 0629-0. 0896-0. 1229-0. 1573-0. 1905-0. 2 279-0. 26 13-0. 2883-0. 3 171-0. 3373-0. 34 77-0. 3584 

119. 00- 0. 0039-0. 0131-0. 028 8-0. C536-0. 085 1-0. 1212-0. 1663-0. 2 128-0. 2576-0. 3 084-0. 3534-0. 3900-0. 4 289-0. 4563-0. 4704-0. 48 48 

126. 00- 0. 0047-0. 01 59-0. 0349-0. 0650-0. 1032-0. 1470-0. 2016-0. 2580-0. 31 24-D. 3739-0. 42B5-0. 4728-0. 5200-0. 5532-0. 5703-0. 5878 

133.00- 0.0055-0.0184-0.0405-0.0754-0.1198-0.1705-0.2339-0.2994-0.3624-0.4338-0.4972-0.5486-0.6034-0.6419-0.6617-0.6820 

139.00- 0.0060-0.0204-0.0448-0.0835-0.1325-0.1887-0.2589-0.3313-0.4011-0.4800-0.5502-0.6071-0.6677-0.7103-0.7323-0.7547 

145.00- 0.0066-0.0221-0.0486-0.0906-0.1439-0.2048-0.2810-0.3596-0.4353-0.5210-0.5972-0.6589-0.7247-0.7710-0.7948-0.8192 

150. 00- 0. 0069-0. 0234-0. 0514-0. C958-0. 152 1-0. 2165-0. 2970-0. 3802-0. 4602-0. 5508-0. 631 3-0. 6966-0. 7662-0. 815 1-0. 8403-0. 8660 

155.00- 0.0073-0.0245-0.0538-0. 1002-0 . 1592-0. 2266-0 . 3 109-0. 3978-0 . 4816-0. 5 76 5-0. 6607-0. 7290-0. 8018-0. 8530-0 . 8794-0 . 9063 

160.00- 0.0075-0.0254-0.0557-0. 1039-0. 1650-0.2350-0.3223-0.4125-0.4994-0.5977-0.6850-0.7558-0.8314-0.8844-0.9118-0.9397 

164.00- 0.0077-0.0260-0.0570-0. 1063-0. 1688-0 . 2404-0 . 3297-0 . 4220-0 . 5109-0 . 6 114-0. 7008-0. 7732-0. 8 505 -0. 9047-0 .93 27-0 . 9613 

168.00- 0.0078-0.0264-0.0580-0. 1082-0. 1718-0.2446-0. 3355-0.4294-0. 5198-D .6222-0.7131-0.7868-0.8654-0.9206-0. 9491-0.9781 

171.00- 0.0079-0.0267-0.0586-0. 1092-0. 1 735-0. 2470-0 .3388-0. 4336-0. 5249-0. 6282-0. 7200-0. 7945-0. 8738-0. 9296-0 .9 584-0. 9877 

174.00- 0.0080-0.0269-0.0590-0. 1100-0 . 1747-0. 2487-0. 341 1-0. 4366-0. 5285- 0.6326-0. 7250-0. 8000-0. 8799-0. 9360-0 .9650-0. 9945 

176.00- 0.008 0-0.0269-0.0592-0. 1103-0. 1752-C. 2494-0. 3422-0. 4379-0. 530 1-0. 6345-0. 7272-0. 8024-0. 8 826-0. 9389-0. 9 679-0. 9976 

177.00- 0.0080-0.0270-0.0592-0.1104-0.1754-0.2497-0.3425-0.4384-0.5307-0.6352-0.7280-0.8033-0.8835-0.9399-0.9690-0.9986 

178.00- 0.0080-0.0270-0.0593-0. 1105-0. 1755-0.2499-0.3428-0.4387-0. 5311-0.6357-0.72 86-0. 8039-0.8842-0.9406-0 .96 S7-0. 99 94 

179. 00- 0. 0080-0. 02 70-0. 0593-0. 1106-0. 1756-0. 2500-0. 3429-0. 4389-0. 5314-0. 6360-0. 7289-0. 8042-0. 8846-0. 9410-0 .9702-0. 9998 

Co 180.00-0.008D-0.0270-0.0593-0. 1106-0. 1756-0. 25C0-0. 3430-0. 4390-0. 5314-0.6361-0. 7290-0.8044-0.8847-0.9412-0.9703-1.0000 

W XCTR= 0. 



00 


RESIDUALS- RHO 

THETA 0.20D 0.300 0.390 0.480 0.560 0.630 0.700 0.760 0.810 0.860 0.900 0.930 0.960 0.980 0.990 1.000 

45.00 0. 0. 0* 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0* 0. 

56.00- 0.0001 0.0000 0 . 0000-0. 0000-0 . 000 1-0. 000 1-0 . 000 l-o. 0002-0 . 0002-0. 0002-0. 000 1-0. 0001 -0. 0001-0. 0000-0 .00 CO 0. 

66. 00- 0. 0032-0. 0001-0. 0000-0. 000 1-0. 0001-0. 0001-0. 0002-0. 0002-0. 000 1-0. 300 1-0. 000 1-0. 0001-0. 0001 -0.0000-0.0 i 0 0. 

76. 00- 0. OOO 3-0. 0002-0. 0002-0. 0002-0. 0002-0. 0002-0. 0002-0. 0002-0. 0002-0. 000 1-0. 000 1-0. 000 1-0. 0000-0. 0000-0. 00 00 0. 

86. 00- 0. 0004-0. 0004-0. 0004-0. 0003-0. 0003-0. 0003-0. 0003-0. 0003-0. 0002- 0. 3002-0. 0001-0. 000 1-0. 0 001-0. 0000-0. 00 00 0. 

95. 00- 0. 0006-0. 0006-0. 000 5-0. 0005-0. 0004-0. 0003-0. 0003-0. 0002-0. 0002-0. 0002-0. 000 1-0. 0001 -0.0001 -0.0000-0.0 00-0. 

103. 00- 0. 003 7-0. 0007-0. 000 7-0. 0005-0. 0004-0. 0004-0. 0003-0. 0002-0. 0001 -0. 0001-0. 0001-0. 0000-0. 0000 -0.0000-0. OOCO-O. 

1 11. 00- 0. 000 B-0. 0009-0. 000 8-0. C006-0. 0005-0. 0004-0. 0002-0. 0002-0. 000 1-0. 0001-0. 0000-0. 0000-0. 0000 0.0000 0.0000-0. 

119. 00- 0. 001 0-0. 00 11-0. 0010-0. 0007-0. 0005-0. 0004-0. 0003-0. 0002-0. 0001-0. 0001-0. 0000-0. 0000-0. 0000-0. 0000-0 .OOCO-O. 

126. 00- 0. 0011-0. 0012-0. 0011-0. C008-0. 0006-0. 0004-0. 0002-0. 0001-0. 0001-0. 0000 0.0000 0.0000 0.0000 0.0000 O.OGCu-O. 

133.00- 0.0011-0.0013-0.0012-0.0009-0.0006-0.0005-0.0002-0.0001-0.0001-3.0000 0.0000 0.0000 0.0000 0.0000 0. OOCO-O. 

139.00- 0.0012-0.0014-0.0013-0.0009-0.0006-0.0005-0.0002-0.0001-0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000-0. 

145.00- 0.0013-0.0015-0.0013-0.0009-0.0006-0.0005-0.0002-0.0001-0.0000 0.0000 0.0001 0.0000 0.0000 0.0000 0.0 CO-O. 

150.00- 0.0013-0.0015-0.0014-0.0010-0.0006-0.0005-0.0002-0.0001-0.0003 0.3001 0.0001 0.0001 0.0001 0.0000 0. OOCO-O. 

155.00- 0.0014-0.0016-0.0014-0.0010-0.0007-0.0005-0.0002-0.0000 0.0000 0.0001 0.0031 0.0001 0.0001 0.0000 0.0000-0. 

160.00- 0.0014-0.0016-0.0015-0.0010-0.0007-0.0005-0.0002-0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0. 001,0-0. 

164. 00- 0. 0014-0. 0016-0. 0015-0. 0010-0. 0007-0. 0005-0. 0002-0. 0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0.0000-0. 

168.00- 0.0014-0.0016-0.0015-0.0010-0.0007-0.0005-0.0002-0.0000 0.0000 0.3001 0.0001 0.0001 0.0001 0.0000 0. OOCO-O. 

171. 00- 0. 0014-0. 0017-0. 0015-0. C010-0. 0007-0. 0005-0. 0002-0. 0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0. OOCO-O. 

174. 00- 0. 0014-0. 0017-0. 0015-C. 0010-0.0007-0.0005-0.0002-0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0.0000-0. 

176. 00- 0. 0014-0. 0017-0. 0015-0. 0010-0. 0007-0. 0005-0. 0002-0. 0000 0.0001 0.0001 0.0001 0.0001 0.0001 0.0000 O.OOCO-O. 

177. 00- 0. 0014-0. 0017-0. 0015-0. 0010-0. 0007-0. 0005-0. 0002 0.0000 0.0001 0.3001 0.0001 0.0001 0.0001 0.0000 O.O.cO-O. 

178. 00- 0. 0014-0. 00 17-0. 001 5—0. 30 10-0.000 7- 0.0005 -0.0002 0.0000 0.0001 0.0001 0.0001 0.0001 0.0001 0.0000 0.000O-0. 

179.00- 0.3015-0.0017-0.0015-0.0010-0.0007-0.0005-0.0002 0.0000 0.0001 0.0001 0.0001 0.0001 0.0001 0.0000 0.0000-0. 

180.00- 0.0015-0.0017-0.0015-0.0010-0.0007-0.0005-0.0002 0.0000 0.0001 0.3001 0.0001 0.0001 0.0001 0.0000 O.O'bO-O. 



TEST CASE 3, THET A0 = 60 DEGREES. 
(WITH AXIAL BOUNDARY VALUES). 

AXIAL VALUES- 
R V 

0 . 0 . 

0.200 -0.008 
0.300 -0.027 

0.390 -0.059 

0.480 -0.111 

0.560 -0.176 

0.630 -0.250 

0.700 -0.343 

0.760 -0.439 

0.810 -0.531 

0.860 -0.636 

0.900 -0.729 

0.930 -0.804 

0.960 -0.885 

0.980 -0.941 

0.990 -0.970 

1.000 - 1.000 


INPUT- 

THETA 

0.200 

0.300 

0.390 

0.480 

0.560 

0.630 

0.700 

RHO 

0.760 

0.810 

0.860 

0.900 

0.930 

0.960 

0.980 

C.990 

l.OJO 

60.00 

0.00 

0.01 

0.03 

0.06 

0.09 

0.13 

0.17 

0.22 

0.27 

0.32 

0.36 

0.40 

0.44 

0.47 

0.49 

0.50 

70.00 

0.68 

1.03 

1.33 

1.64 

1.92 

2.15 

2.39 

2.60 

2.77 

2.94 

3.08 

3. 18 

3.28 

3.35 

3.39 

0.34 

79.00 

0.38 

0.57 

0.74 

0.92 

1.07 

1.20 

1.34 

1.45 

1.55 

1.64 

1.72 

1.77 

1.83 

1.87 

1.69 

0.19 

88.00 

0.07 

0.10 

0.14 

U.17 

0.20 

0.22 

0.24 

0.27 

0.28 

0.30 

0.31 

0.32 

0.34 

0.34 

0.35 

0.03 

96.00 

-0.21 

-0.31 

-0.41 

-0.50 

-0.59 

-0.66 

-0.73 

-0.79 

-0.85 

-0.90 

-0.94 

-0.97 

-1.00 

-1.02 

-1.03 

-0.10 

104.00 

-0.48 

-0.73 

-0.94 

-1.16 

-1.35 

-1.52 

-1.69 

-1.84 

-1.96 

-2.08 

-2.18 

-2.25 

-2.32 

-2.37 

-2.40 

-0.24 

112.00 

-0.75 

-1.12 

-1.46 

-1.80 

-2.10 

-2.36 

-2.62 

-2.85 

-3.03 

-3.22 

-3.37 

-3.48 

-3.60 

-3.67 

-3.71 

-0.37 

119.00 

-0.97 

-1.45 

-1.89 

-2.33 

-2.71 

-3.05 

-3.39 

-3.68 

-3.93 

-4.17 

-4.36 

-4.51 

-4.65 

-4.75 

-4. 80 

-0.48 

126.00 

-1.18 

-1.76 

-2.29 

-2.82 

-3.29 

-3.70 

-4.11 

-4.47 

-4.76 

-5.05 

-5.29 

-5.47 

-5.64 

-5.76 

-5.82 

-0.59 

132.00 

-1.34 

-2.01 

-2.61 

-3.21 

-3.75 

-4.22 

-4.68 

-5.09 

-5.42 

-5.75 

-6.02 

-6.22 

-6.42 

-6.56 

-6.62 

-0.67 

138.00 

-1.49 

-2.23 

-2.90 

-3.57 

-4.16 

-4.68 

-5.20 

-5.65 

-6.02 

-6.39 

-6.69 

-6.91 

-7.13 

-7.28 

-7.36 

-0.74 

144.00 

-1.62 

-2.43 

-3.16 

-3.88 

-4.53 

-5.10 

-5.66 

-6.15 

-6.55 

-6.96 

-7.28 

-7.52 

-7.77 

-7.93 

-8.01 

-0. 81 

149.00 

-1.71 

-2.57 

-3.34 

-4.11 

-4.80 

-5.40 

-6.00 

-6.51 

-6.94 

-7.37 

-7.71 

-7.97 

-8.23 

-8.40 

-8.49 

-0.86 

154.00 

-l.BO 

-2.70 

-3.51 

-4. 31 

-5.03 

-5.66 

-6.29 

-6.83 

-7.28 

-7.73 

-8.09 

-8.36 

-8.63 

-8.81 

-8. 90 

-C.90 

158.00 

-1.85 

-2.78 

-3.62 

-4.45 

-5.19 

-5.84 

-6.49 

-7.05 

-7.51 

-7.97 

-8.34 

-8.62 

-8.90 

-9.09 

-9.18 

-0.93 

162.00 

-1.90 

-2.85 

-3.71 

-4.57 

-5.33 

-5.99 

-6.66 

-7.23 

-7.70 

-8.18 

-8.56 

-8.84 

-9.13 

-9.32 

-9. 42 

-0.95 

166.00 

-1.94 

-2.91 

-3.78 

-4.66 

-5.43 

-6.11 

-6.79 

-7.37 

-7.86 

-8.34 

-Q.73 

-9.02 

-9.31 

-9.51 

-9.61 

-0.97 

169.00 

-1.96 

-2.94 

-3.83 

-4.71 

-5.50 

-6.18 

-6.87 

-7.46 

-7.95 

-8.44 

-8.83 

-9.13 

-9.42 

-9.62 

-9.72 

-0.98 

172.00 

-1.98 

-2.97 

-3.86 

-4.75 

-5.55 

-6.24 

-6.93 

-7.53 

-8.02 

-8.52 

-8.91 

-9.21 

-9.51 

-9.70 

-9. 80 

-0. 99 

174.00 

-1.99 

-2.98 

-3.88 

-4.77 

-5.57 

-6.27 

-6.96 

-7.56 

-8.06 

-3.55 

-8.95 

-9.25 

-9.55 

-9.75 

-9.85 

-0.99 

176.00 

-2.00 

-2.99 

-3.89 

-4.79 

-5.59 

-6.28 

-6.98 

-7.58 

-8.08 

-8.58 

-8.98 

-9.28 

-9.58 

-9.78 

-9.88 

-1.00 

177.00 

-2.00 

-3.00 

-3.89 

-4.79 

-5.59 

-6.29 

-6.99 

-7.59 

-8.09 

-8.59 

-8.99 

-9.29 

-9.59 

-9.79 

-9.89 

-1.00 

178.00 

-2.00 

-3.00 

-3.90 

-4.80 

-5.60 

-6.30 

-7.00 

-7.60 

-8.10 

-8.59 

-8.99 

-9.29 

-9.59 

-9.79 

-9.89 

-1.00 

179.00 

-2.03 

-3.00 

-3.90 

-4.80 

-5.60 

-6.30 

-7.00 

-7.60 

-8.10 

-8.60 

-9.00 

-9.30 

-9.60 

-9.80 

-9.90 

-1.00 

180.00 

-0.01 

-0.03 

-0.06 

-0.11 

-0.18 

-0.25 

-0.34 

-0.44 

-0.53 

-0.64 

-0.73 

-0.80 

-0.88 

-0.94 

-0. 97 

-1.00 


00 

on 



w OUTPUT- RHO 

CT> THETA 0.203 0.300 0.390 0.480 0.560 0.630 0.700 0.760 0.810 0.860 0.900 0.930 0.960 0.980 C.990 1.000 

60.00 0.0040 0.0135 0.0297 0.u553 0.0878 0.1250 0.1715 0.2195 0.2657 0.3180 0.3645 0.4022 0.4424 0.4706 0.4851 0.5000 

70.00 0.0026 0.0091 0.0202 0.0377 0.0600 0.0854 0.1172 0.1500 0.1816 0.2174 0.2492 0.2750 0.3025 0.3219 0.3318 0.3420 

79.00 0.0013 0.0049 0.0111 0.0209 0.0333 0.0475 0.0653 0.0836 0.1013 0.1212 0.1390 0.1534 0.1688 0.1796 0.1851 0.1908 

88.00- 0.0001 0.0006 0.0017 0.0036 0.0059 0.0085 0.0118 0.0151 0.01B4 0.0221 0.0253 0.0280 0.0308 0.0328 0.0338 0.0349 

96.00- 0. 001 3-0. 0033-0. 006 7-0. 01 19-0. 0187-0. 0264-0. 0360-0. 0460-0. 0557-0. 0666-0. 0763-0. 0841-0. 0925-0. 0984-0. 10 14-0. 1045 

104.00- 0.0025-0.0072-0.0150-0.0272-0.0428-0.0608-0.0832-0. 1063-0.1287-0.1539-0.1764-0. 1946-0.2140-0.2277-0.2347-0.2419 

1 12. 00- 0. 0037-0. 01 09-C. 0230-0. C420-0. 0662-0. 0940-0. 128 7-0. 1646-0. 1992-0. 2 383-0. 2731-0. 30 13-0. 3 314-0. 3526-0. 3635-0. 3746 

119.00- 0.0047-0.0140-0.0296-0.0542-0.0856-0.1216-0.1665-0.2129-0.2577-0.3084-0.3534-0.3900-0.4289-0.4563-0.4704-0.4848 

126.00- 0.0056-0.0169-0.0358-0.0657-0.1037-0.1473-0.2018-0.2581-0.3124-0.3739-0.4285-0.4728-0.5200-0.5532-0.5703-0.5878 

132.00- 0.0063-0.0192-0.0407-0.0747-0. 1180-0.1677-0.2297-0.2938-0.3556-0.4256-0.4877-0.5382-0.5920-0.6298-0.6492-0.6691 

138.00- 0.0069-0.0212-0.0452-0.0829-0. 1310-0.1862-0.2551-0.3263-0.3949-0.4726-0.5417-0.5977-0.6574-0.6994-0.7211-0.7431 

144. 00- 0. 0075-0. 02 31-0. 049 1-0. 0903-0. 142 6-0. 2027-0. 2776-C. 3552-0. 4299-0. 5 145-0. 5897-0. 6507-0. 7 157-0. 7614-0. 78 50-0. 8090 

149.00- 0.0079-0.0244-0.0520-0.0956-0.1511-0.2147-0.2941-0.3763-0.4555-0.5451-0.6248-0.6894-0.7583-0.8067-0.8317-0.8572 

154.00- 0.0082-0.0255-0.0545-0.1002-0.1584-0.2251-0.3084-0.3945-0.4776-0.5716-0.6551-0.7229-0.7951-0.8459-0.8721-0.8988 

158.00- 0.0085-0.0263-0.0562-0. 1033-0.1633-0.2322-0.3181-0.4070-0.4927-0.5896-0.6758-0.7457-0.8202-0.8726-0.8996-0.9272 

162.00- 0.0086-0.0,269-0.0576-0. 1060-0.1675-0.2381-0.3263-0.4174-0.5054-0.6048-0.6932-0.7649-0.8414-0.8951-0.9228-0.9511 

166.00- 0.0088-0.0274-0.0587-0. 1081-0.1709-0.2429-0.3329-0.4259-0.5156-0.6170-0.7072-0.7804-0.8584-0.9132-0.9415-0.9703 

169.00- 0.0088-0.0277-0.0593-0. 1093-0 . 1 728-0. 2458-0 . 3368-0. 4309-0 . 5216-0 . 6242-0. 7155-0 . 7895-0 . 8684-0. 9239-0 .95 24-0 . 9816 

172.00- 0. JOB 8-0. 02 78-0. 0598-0. 1102-0.1743-0.2479-0.3397-0.4346-0.5262-0.6297-0.7218-0.7964-0.8 760-0.9320-0.9608-0.9903 

174.00- 0.0088-0.0279-0.0600-0. 1106-0 . 1750-0 . 2489-0. 341 2-0. 4365-0. 52 84- 0 . 6 324-0. 7249-0. 7998-0. 8798-0. 9360-0 .965 U-0 . 9945 

176.00- 0.0087-0.0279-0.0600-0. 1109-0.1755-0.2497-0.3422-0.4378-0.5301-0.6344-0.7271-0.8023-0.8825-0.9389-0.9679-0.9976 

177.00- 0.0087-0.0278-0.0600-0. 11 10-0 . 1757-0. 2499-0. 3426-0. 4383-0 . 5306-0 . 6351-0. 7279-0 . 8032-0 . 8835-0. 9399-0 .96 89-0 .9986 

178.00- 0.0086-0.0277-0.0600-0. 1110-0.1758-0.2501-0.3428-0.4387-0.5310-0.6356-0.7284-0.8038-0.8841-0.9406-0.9697-0.9994 

179.00- 0.0084-0.0276-0.0598-0. 1109-0 . 1 758-0. 2502-0 . 3430-0. 4389-0. 5313- 0 . 6 359-0. 7288-0. 8042-0. 8 846-0. 9410-0 .9 701-0 . 9998 

180.00- 0.0080-0.0270-0.0593-0. 1106-0.1756-0.2500-0.3430-0.4390-0.5314-0.6361-0.7290-0.8044-0.8847-0.9412-0.9103-1.0000 
XCTR* 0. 

EXACT ANSWERS RHO 

THETA 0.200 0.300 0.390 0.480 0.560 0.630 0.700 0.760 0.810 0.860 0.900 0.930 0.960 0.980 C.990 1.000 

60.00 0.0040 0.0135 0.0297 0.0553 0.0878 0.1250 0.1715 0.2195 0.2657 0.3180 0.3645 0.4022 0.4424 0.4706 0.4851 0.5000 

70.00 0.0027 0.0092 0.0203 0.0378 0.0601 0.0855 0.1173 0.1501 0.1818 0.2175 0.2493 0.2751 0.3026 0.3219 0.3319 0.3420 

79.00 0.0015 0.0052 0.0113 0.C211 0.0335 0.0477 0.0654 0.0838 0.1014 0.1214 0.1391 0.1535 0.1688 0.1796 0.1851 0.1908 

88.00 0.0003 0.0009 0.0021 0.0039 0.0061 0.0087 0.0120 0.0153 0.0185 0.0222 0.0254 0.0281 0.0309 0.0328 0.0i39 0.0349 

96.00- 0.0008-0.0028-0.0062-0.0116-0.0184-0.0261-0.0359-0.0459-0.0556-0.0665-0.0762-0.0841-0.0925-0.0984-0.10 14-0. 1045 

104. 00- 0. 001 9-0. 0065-0. 0144-0. 0268-0. 0425-0. 0605-0. 0830-0. 1062-0. 1286-0. 1 539-0. 1764-0. 1946-0. 2 140 -0.2277-0. 2347-0. 24 19 

112.00- 0.0030-0.0101-0.0222-0.0414-0.0658-0.0937-0.1285-0.1644-0.1991-0.2383-0.2731-0.3013-0.3314-0.3526-0.3635-0.3746 

119. 00- 0. 003 9-0. 01 31-0. 028 8-0. C536-0. 085 1-0. 1212-0. 1663-0. 2128-0. 2 576- 0.3084-0. 3534-0. 3900-0. 4289-0. 456 3-0 .4704-0. 48 48 

126.00- 0.0047-0.0159-0.0349-0.0650-0.1032-0.1470-0.2016-0.2580-0.3124-0.3739-0.4285-0.4728-0.5200-0.5532-0.57:3-0.5878 

132.00- 0.0054-0.0181-0.0397-0.0740-0.1175-0.1673-0.2295-0.2937-0.3556-0.4256-0.4878-0.5382-0.5920-0.6298-0.6493-0.6691 

138.00- 0. 9059-0.0201-0.0441-0.0822-0.1305-0.1858-0.2549-0.3262-0.3949-0.4727-0.5418-0.5978-0.6575-0.6994-0.7211-0.7431 

144.00- 0.0065-0.0218-0.0480-0.0895-0.1421-0.2023-0.2775-0.3551-0.4299-0.5146-0.5898-0.6507-0.7158-0.7614-0.7850-0.8090 

149.00- 0.0069-0.0231-0.0508-0.0948-0.1505-0.2143-0.2940-0. 3763-0. 4555-0. 545?-0. 6249-0. 6895-0.7584-0.8068-0.8317-0.8572 

154. 00- 0. 0072-0. 0243-0. 053 3-0. 0994-0. 1578-0. 2247-0. 3083-0. 3945-0. 4777-0. 5 717-0. 6552-0. 72 30-0. 7 952-0. 8459-0. 87 21-0. 8988 

158.00- 0.0074-0.0250-0.0550-0. 1025-0.1628-0.2318-0.3180-0.4070-0.4927-0.5897-0.6759-0.7458-0.8203-0.8727-0.8996-0.9272 

162.00- 0.0076-0.0257-0.0564-0. 1052-0.1670-0.2378-0.3262-0.4175-0.5054-0.5049-0.6933-0.7650-0.8414-0.8951-0.9228-0.9511 

166.00- 0.0079-0.0262-0.0576-0. 1073-0 . 1704-0. 2426-0. 3328-0 . 4259-0. 51 57-0 . 5 172-0. 7073-0. 7805-0. 8 585-0. 9132-C .94 15-0 . 9703 

169.00- 0.0079-0.0265-0.0582-0. 1086-0.1724-0.2455-0.3367-0.4309-0.5217-0.6244-0.7156-0.7896-0.8685-0.9239-0.9525-0.9816 

172.00- 0.0079-0.0267-0.0587-0. 1095-0.1739-0.2476-0.3397-0.4347-0. 5263-0.6299-0. 7219-0. 7965-0. 8 761-0. 9320-0.96C9-0. 9903 

174.00- 0.0080-0.0269-0.0590-0.1100-0.1747-0.2487-0.3411-0.4366-0.5285-0.6326-0.7250-0.8000-0.8799-0.9360-0.9650-0.9945 

176.00- 0.0080-0.0269-0.0592-0. 1103-0. 1752-0.2494-0. 3422-0. 4379-0. 5301-0.6345-0.7272-0. 8024-0. 8826-0.9389-0.9679-0.9976 

177.00- 0.0080-0.0270-0.0592-0. 1104-0.1754-0.2497-0.3425-0.4384-0.5307-0.6352-0.7280-0.8033-0.8835-0.9399-0.9690-0.9986 

178.00- 0.0080-0.0270-0.0593-0. 1105-0.1755-0.2499-0.3428-0.4387-0.5311-0.6357-0.7286-0.8039-0.8842-0.9406-0.9697-0.9994 

179.00- 0. 0080-0.0270-0.059 3-0. 1106-0. 1756-0. 2500-0. 3429-0.4389-0. 5314- 0.6360-0. 7289-0. 8042-0. 8 846-0. 9410-0 .9702-0.99 98 

180. 00- 0. 0080-0. 0270-0. 0593-0. 1106-0. 1756-0. 2500-0. 3430-0. 4390-0. 5314-0. 6361-0. 7290-0. 8044-0. 8 847-0. 9412-0. 9703-1. 0000 
XCTR* 0. 



RESIDUALS- RHO 

THETA 0.233 0.300 0.390 0.480 0.560 0.630 0.700 0.760 0.813 0.860 0.900 0.930 0.960 0.980 0.990 1*000 

60.00 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0* 0. 0. 

70. 00- 0. 003 1-0. 0001-0. 000 1-0. 000 1-0. 000 1-0. 000 1-0. 000 1-0. 0001-0. 000 1-0. 300 1-0. 000 1-0. 00 01-0* 0001-0. 0000-0. 00 CO 0. 

79. 00- 0. 0002 -0.0002-0. 000 2-0. 0002 -0.0002-0. 0002-0. 0002-0. 0002-0. 000 1-3. 0001-0. 000 1-0. 0001-0. 0000-0. 0000-0. 00 0. 

88. 00- 0. 0034-0. 0004-0. 0003-0. 0003-0. 0003-0. 0002-0. 0002-0. 0002-0. 0002-0. 0001-0. 0001-0. 0001-0. 0000 -0.0000-0.0 -‘CO 0. 

96. 00- 0. 0005-0. 0005-0. 000 5-0. 0004-0. 0003-0. 0003-0. 0002-0. 0001-0. 000 1-0. 000 1-0. 0001-0. 0000-0. 0 000 -0.0000-0. 00 CO-O. 

104. 00- 0. 0006-0. 0007-0. 0006-0. 000 5-0. 0004-0. 0003-0. 0002-0. 0001-0. 0001-0. 0001-0. 0000-0. 00 00-0. 0000-0. 0000-0. 00 GC-O. 

1 12. 00- 0. 0007-0. 0008-0. 0008-0. 0006-0. 0004-0. 0003-0. 0002-0. 0001-0. 000 1-0. 000 1-0. 0000-0. 0000-0. 0000-0. 0000- G. 0000-0. 

119. 00- 0. 0008-0. 0009-0. 0009-0. 0006-0. 0005-0. 0004-0. 0002-0. 0001-0. 0001-0. 0000 0.0000 0.0000 0.0000 0.0000 O.OOuO-O. 

126.00- 0.0009-0.0010-0.0010-0.0007-0.0005-0.0004-0.0002-0.0001-0.0001-0.0000 0.0000 0.0000 0.0000 0.0000 0.0300-0. 

132. 00- 0. 0009-0. 0011-0. 0010-0. U007-0. 0005-0. 0004-0. 0002-0. 0001-0. 0000 0.0000 0.0001 0.0000 O.OCOO 0.0000 0 .0000-0. 

138.00- 0.0010-0.0012-0.0011-0.0008-0.0005-0.0004-0.0002-0.0000-0.0000 0.0001 0.0001 0.0001 0.0000 0.0000 C. 000^-0. 

144. 00- 0. 0010-0. 0012-0. 0011-0. 0008-0. 0005-0. 0004-0. 0001-0. 0000 0.0000 0.0001 O.OOGl 0.0001 0.0000 0.0000 O.OOCC-O. 

149.00- 0.0010-0.0013-0.0012-0.0008-0.0005-0.0004-0.0001-0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 C.OOCu-O. 

154.00- 0.0010-0.0013-0.0012-0.0008-0.0005-0.0004-0.0001 0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0.0 GC-O. 

158. 00- 0. 0013-0. 0013-0. 0012-0. G008-0. 0005-0. 0003-0. 0001 0.0000 0.0001 3.0001 0.0031 0.0001 0.0001 0.0000 C.OO^C-O. 

162.00- 0.0010-0.0013-0.0012-0.0008-0.0005-0.0003-0.0001 0.0000 0.0001 0.0001 0.0001 0.0001 0.0001 0.0000 Q.OOOu-O. 

166.00- 0.0010-0.0012-0.0011-0.0008-0.0005-0.0003-0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0000 0.0 cO-0. 

169. 00- 0. 0039-0. 0012-0. 0011-0. 0007-0. 0005-0. 0003-0. GuOl 0.0001 0.0001 0.3001 0.0001 0.0301 0.0001 0.0000 O.OQOj-O. 

172.00- 0.0039-0.0011-0.0010-0.0007-0.0004-0.0003-0.0000 0.0001 0.0031 0.3001 0.0001 0.0001 0.0001 0.0000 O.OOCC-O. 

174.00- 0.0008-0.0010-0.0010-0.0006-0.0004-0.0003-0.0000 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0000 0.0000-0 • 

176.00- 0.0037-0.0009-0.0009-0.0006-0.0004-0.0002-0.0000 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0000 O.GQQc-O. 

177.00- 0.0007-0.0008-0.0008-0.0005-0.0003-0.0002-0.0000 0*0001 0.0001 0.3001 0.0001 0.0001 0.0001 0.0000 O.OOGC-O. 

178. 00- 0. 0036-0.0007-0.0007-0.0005-0. 0003-0. 0002-0. 0000 0.0001 0.0001 0.0001 0.0031 0.0001 0.0001 0.0000 0.0000-0. 

179.00- 0.0004-0.0006-0.0005-0.0003-0.0002-0.0001-0.0000 0.0000 0.0001 0.0001 0.0001 0.0001 0.0000 0.0000 O.OGGC-O. 

180.00- 0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. 

*01* EXIT IN TEST 
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APPENDIX B 


LISTING OF FORTRAN PROGRAM 


The following computer program should be considered a subroutine, or black box 
for which an executive program is written. Such a program and sample test run are 
given in appendix A. 

Subroutine POISS1 generates the r-0 mesh then generates and inverts the requisite 
coefficient matrices. POISS2 then uses these matrices together with the right-hand side 
values to generate the solution. 

SIBFTC P0IS1 

SUBROUTINE P0ISS1 
C 

COMMON/ BLOCK 1/NR, NT, R A (17) • TA < 2 5 ) * X( 25* 16) , XI, AVS 
COMMON/BLOCK2/V(25,25,15),A(25,15),FM(25,15),BCMC(25),BCMR(16), 

1 BCMRI (25) , BCMA( 16)* V 1 , F Ml , NRM 1 , NRM2 , NRM3, NTM1 , RBA ( 16 I ,DRA(16> , 

2 NTE*NTEM1,TBA(25)*0TA(25) , S TBA ( 2 5 ) , R OW ( 25 ) ,0(25,25) 

LOGICAL AVS 

C 

NRM 1=NR-1 
NRM2=NR-2 
NR M3=NR-3 
NTM1=NT-1 
C 

C CHECK MESH SIZE. 

IFINR.GT.17) CALL AR ERR ( 23HT 00 MANY R MESH LINES $1 
IF(NT.GT.25) CALL AR ERR ( 27HT00 MANY THETA MESH LINES $) 

C 

C GENERATE RBAR AND DEL R ARRAYS. 

DO 1 I R= 1 , NRM1 

RBA ( I R ) = .5*(RA«IR)fRAl IR+1) ) 

1 DRA ( I R ) = RA( IR + 1 )-RA( IR) 

C 

C GENERATE THETA-BAR, SIN THETA-BAR 

C AND DEL-THETA ARRAYS. 

DO 2 IT=1,NTM1 

TBA (IT)- . 5*(TA(IT)+TA(IT+1) ) 

STB A ( I T ) = S I N ( TB A ( IT) ) 

2 DTA ( I T ) = TA{ IT+1 )-TA ( IT) 

C 

C TEST CONFIGURATION FOR FIRST PASS. 

IF ( TA( 1 ) • EO. 0 • ) GO TO 3 
NTE=NT 

IF(AVS) NT E=NT-1 
NTEM1=NTE-1 
CALL SEG1 
RETURN 
C 

C FULL SPHERE SOLUTION (1ST PASS). 

3 CALL SPHER1 
RETURN 

C 

ENTRY P0ISS2 
IF(TA(l)„EQ.O. ) GO TO 4 
CALL SEG2 
RETURN 
C 

C FULL SPHERE SOLUTION (2ND PASS). 

4 CALL SPHER2 
C 

RETURN 

END 
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SIBFTC SPHR1 


SUBROUTINE SPHER1 

GENERATES AND INVERTS COEFFICIENT MATRICES 
USED FOR FULL SPHERE SOLUTION. 

NR— 2 TRIDIAGONAL BLOCK MATRICES. 

2 ( NR— 2 ) -1 DIAGONAL BLOCK MATRICES. 

COMMON/BLOCK 1/NR ,NT,RA<17) , TA ( 25 ) » X ( 25 , 16 ) f Xl, AVS 

COM MON/ BLOC K2/V( 25*25*15 )*A(25«15)*FM(25«15) #BCMC( 25 ) . BCMR( 16) 

1 BCMR1 (25) $ BCMA ( 16 ) , VI ,FM1 ,NRM1 .NRM2.NRM3.NTM1 ,RBA { 16) f DRA( 16) 

2 NTE'fNTEMl f TBA(25)»DTA(25) « S TBA ( 25 ) . ROW( 25) ,0(25 .25) 

LOGICAL AVS 

INVERSE OF CENTER POINT. 

Vl= 1. / AE ( 1 v ItO) 

FM1= FE ( 1 t 1 ) 

Kl=l 


DO 14 I R a 2 . NRM1 
GO TO ( 1,3),K1 

I R*2 ONLY. 

1 DO 2 I T I s 1 . NT 

A( I T1 t 1 ) = AE(2,IT1,1) 

FM ( IT1 . 1 ) = FEU.ITl) 

2 BCMC(IT1)=-AE(NRM1,IT1,2) 

GO TO 5 

3 DO 4 I T 1= 1 » NT 
AUTl.IR-l)* AE( IRflTlfl ) 

4 FM ( IT1 » IR— 1 ) = FE( IRt IT1) 

GENERATE NR-2 NTXNT MATRICES. 

(ONE FOR EACH R-CONSTANT LINE). 

5 DO 12 I T 1 = I* NT 
GO TOt 6.8) * K1 

FORM -A2*V1*A2. 

6 DO 7 I T2= 1 ♦ NT 

7 D( iTlt IT2)=— A( ITltl) *V1*A( IT2 # 1) 

GO TO 10 

FORM — A( I )*V( I-1)*A(I ). 

8 DO 9 IT2=1,NT 

9 DC I T1 v I T2 ) AC IT1, IR-l)*Vt I T 1, I T2, IR-2 ) * A( IT2.IR-1) 

10 IF(ITi.EQ.l) GO TO 11 

ADD TRIDIAGONAL ELEMENTS OF 
B ( l ) TO -A<I)*V(I-l)*A(l>. 

D(IT1,IT1-1)= AE( IR. IT1.3)+D( IT1.IT1-1) 

11 D(ITl.ITl)- AE( IR» IT1.0)+D( IT1.IT1) 

IF< IT1.EQ.NT) GO TO 12 

D ( I Tl. I Tl + 1 ) - AE( IR. IT1,4)+D(IT1. IT 1 + 1) 

12 CONTINUE 

CALL F ACTOR ( D.RQW . NT ) 

CALL INVERT(D.R0W,NT.V(1.1. IR-1) ) 

GO TO (13.14), K1 

13 Kl=2 

14 CONTINUE 

RETURN 

END 
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$ I BFTC SPHR2 

SUBROUTINE SPHER2 
C 

COMMON/ BLOCK 1/ NR, NT, RA( 17) ,TA( 25 ) , XI 25, 16) , XI, AVS 
C0MM0N/BL0CK2/VI 25,25, 15 ),A( 25 • 15 ) , FM < 25, 15 ) , BCMC ( 2 5 ) , BCMR ( 16 ) , 

1 BCMRK25) ,BCMA(16), V1,FM1,NRM1,NRM2,NRM3,NTM1,RBA{ 16),DRA(16) , 

2 NTE,NTEM1,TBA(25),DTA(25) , STBAI 25 ) , ROW< 25 ) , D( 25 • 25 ) 

LOGICAL AVS 

C 

DOUBLE PRECISION T,Y(25) 

C CALCULATE Gl. 

X 1= V1*FM1*X1 

CALCULATE G2 (IR=2). 

DO 1 I T=1 , NT 

1 Y { I T ) “ EM (IT,1)*X(IT,1)-A( IT,1)*X1 
DO 3 I T 1= 1 , NT 
T=0 • DO 

DO 2 I T2= 1 , NT 

2 T=T+V(IT1*IT2,1)*Y(IT2) 

3 X( I T 1, 1 ) =T 

CALCULATE G3 TO G ( NRM2) • 

DO 6 I R=3 , NRM2 
DO 4 I T~ 1 , NT 

4 Y ( I T ) = FM(IT,IR-U*XC IT, IR-U-AUT, IR-1)*X( IT, IR-2) 

DO 6 I T 1= 1 ,NT 
T=0 • DO 

DO 5 I T2= 1 , NT 

5 T=T + V(IT1, IT2,IR-1)*Y(IT2> 

6 X ( I T 1 , 1 R— 1 ) = T 

CALCULATE G(NRMl) . 

DO 7 I T=1 , NT 

7 Y ( I T ) = FM ( IT, NRM2 ) *X ( I T ,NRM2 ) +BCMC( IT) *X (IT, NR Ml) -A ( I T ,NRM2 ) * 

1 X ( I T , NRM3 ) 

DO 9 I T 1= 1 ,N T 
T= 0.00 
DO 8 I T2= 1 , NT 

8 T=T+V( IT1,IT2,NRM2)*Y( IT2) 

9 X ( I T1 , NRM2 ) =T 

BACKWARD SUBSTITUTION FOR SOLUTION. 

DO 11 I RC= 1 » NRM3 
I R= NRMl-IRC 
DO 11 I T 1 = 1 , NT 
T= 0. DO 

DO 10 I T2= 1 , NT 

10 T = T +V( I T 1 , I T 2 , I R-l ) * A ( IT 2, I R >*X( IT2, IR ) 

11 X(IT1,IR-1)= X ( I T 1 , 1 R-l ) — T 


C CALCULATE CENTER POINT. 

T = 0 • DO 

DO 12 I T 1= 1 , NT 
12 T=T+A( I T1,1)*X(IT1,1 ) 

Xl= XI -VI *T 
C 

RETURN 

END 
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$ I BFTC SEGM1 


SUBROUTINE SEG1 

BOUNDARY VALUE COEFFICIENTS. 

8CMA -AXIAL VALUES ( THETA=P I 1 • 

BCMC -SURFACE VALUES (CIRCUMFERENCE). 

BCMR -RADIAL LINE VALUES ( THETA=CONST • SURFACE). 
BCMR 1- CENTRAL VALUES. 


COMMON/ BLOCK 1/ NR tNT t R A ( IT) ,TA(25),X(25,16) «X1* AVS 

COMMON/ BLOCK2/V(25,25,15),A(25,15)»FM(25,l5) ,BCMC( 25) f BCMR (16) 

1 BCMR 1(25) »BCMA( 16 ) , V 1 , FM1 , NRM1 , NRM2 , NRM3 , NTM1 , RBA ( 16 ) ,DRA(16) 

2 NTE,NTEM1,TBA(25) f DTA(25) , STBA(25) , ROW (25) ,D(25,25) 

LOGICAL AVS 

DO 1 I R=2 1 NRM1 
BCMR ( I R ) = -AE ( IR.2,3) 

1 IF(AVS) BCMA(IR)=-AE( IR,NTE,4) 

DO 3 I T 1* 1 ? NTE 
INITIALIZE. 

DO 2 I T2= 1 » NT E 

2 D ( I T 1, I T2 ) =0 • 0 

BCMC ( I T 1 ) =-AE ( NR Ml » I T 1 v 2 ) 

BCMR1(IT1)=-AE(2,IT1,1) 

FM( ITlt 1 )-FE( 2, ITl ) 

IF ( IT1 • EQ • 1 ) GO TO 3 

GENERATE ELEMENTS OF 61 THEN INVERT TO GET VI. 

( T, D« B COEFFS. FOR IR=2 MESH LINE). 

IF( IT1.NE.2) D( IT1-1, I T 1-2 ) = AE ( 2 , I Tl , 3 ) 

D( IT1-1, IT1-1)=AE(2, I T 1 * 0 ) 

IF( IT1.NE.NTE) D ( I Tl - 1 , I Tl ) = AE ( 2 • I Tl , 4 > 

3 CONTINUE 

CALL FACTOR(D*ROW»NTEMl) 

CALL INVERT(D,ROW,NTEMl, V( 1, 1, 1) ) 

DO 7 I R=3 # NRM1 
DO 4 I T 1*2 tNTE 
A ( I T1 1 1 R— 1 ) = AE( IR, IT 1 » 1 ) 

4 FM( ITl* IR-1)=FE( IR • I T 1 ) 

GENERATE REMAINING VI MATRICES. 

DO 6 I T 1*2 * NTE 
DO 5 I T2*2» NTE 

5 D( IT1-1, IT2-1)=-A( ITlt I R-l ) * V ( I Tl- 1 , I T2- 1 • I R-2 ) * A ( IT2. IR-1) 

IF( IT1.NE.2) D ( I T 1-1 , IT1-2)=D( I Tl- 1 , I Tl-2 ) +AE ( I R , I T 1 » 3 ) 
D(IT1-1, IT1-1)=D( IT1-1, IT1-1 )+AE( IR, IT1,0) 

6 IF ( ITl. NE. NTE) D ( I Tl- 1 , I Tl ) =D ( I Tl-1 ♦ I T 1 ) +AE ( I R , I Tl , 4 ) 

CALL FACTOR! D, ROW, NTEM1 ) 

7 CALL INVERT(D, ROW, NTE Ml ,V(1,1,I R— 1 ) ) 

RETURN 

END 
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$ I BFTC SEGM2 


SUBROUTINE SEG2 
C 

C0MM0N/BL0CK1/NR,NT,RA(17) , TA(25), X(25, 16), XI, AVS 

COMMON/ BLOCK 2/ V (25*25? 15 )« A( 25* 15) . FM ( 25 , 15) ,BCMC(25) , BCMR ( 16) 

1 BCMR1(25),BCMA(16), V 1 , FMl , NRM1 , NRM2 *NRM3, NTM1 , RBA { 16),DRA(16) 

2 NTE,NTEM1,TBA(25) ,DTA(25) ,STBA(25),R0W(25>.D(25,25> 

LOGICAL AVS 

C 

DOUBLE PRECISION T,Y(25) 

CALCULATE GL (INNERMOST RING). 

DO 1 I T 1 = 2 ♦ NTE 

1 Y( IT1)=FM( IT1, 1)*X(IT1, 1 KBCMRK IT1 )*X1 
Y(2)=Y(2)+BCMR(2)*X( 1,1) 

IF (AVS) Y(NTE)=Y(NTE )+BCMA(2 )*X(NT,1) 

DO 3 I T 1=2 . NTE 
T=O.DO 

DO 2 I T2=2 » NTE 

2 T=T+V ( I Tl-1 , IT 2- 1,1) * Y ( I T2 ) 

3 X ( I Tl. 1 )*T 

CALCULATE G(IR-l) (INTERIOR 
RINGS 2 TO NRM3 ) • 

DO 6 I R=3 . NRM2 
DO 4 I T 1=2 .NTE 

4 Y ( I Tl ) *FM ( I Tl * IR-1)»X( IT It IR-1)-A( I T 1 , I R-l ) *X ( I T 1 , I R-2 ) 
Y(2)=Y(2)+BCMR( IR )*X(1,IR-1) 

IF (AVS) Y(NTE)=Y(NTE ) +BCMA ( I R ) *X ( NT, I R-l) 

DO 6 I T 1=2 .NTE 
T=0 • DD 

DO 5 I T 2=2 , NTE 

5 T=T+V( IT1-1,IT2-1,IR-1)*Y( IT2) 

6 X ( I T 1 , I R— 1 ) = T 

CALCULATE G(NR-2). 

DO 7 IT1=2,NTE 

7 Y( IT1) =FM( IT1 ,NRM2 )*X( IT 1.NRM2 ) +BCMC ( IT1 )*X( IT1.NRM1) 

1 — A ( I T 1 ,NRM2)*X( IT1,NRM3) 

Y( 2)=Y (2)+BCMR(NRMl) *X{ 1.NRM2) 

IF (AVS) Y ( NT E ) =Y ( NTE ) +BCMA ( NRM 1 ) *X ( NT , NRM2 ) 

DO 9 IT1=2,NTE 
T=0 . DO 

DO 8 I T 2= 2 » NTE 

8 T = T + V( IT1-1, IT2-1,NRM2)*Y( IT2) 

9 X( IT1,NRM2)=T 

BACKWARD SUBSTITUTION FOR SOLUTION. 

DO 11 I RC= 1 , NRM3 
I R = NRM1 — I RC 
DO 11 I Tl =2 , NT E 
T=0. DO 

DO 10 I T2 = 2 . NTE 

10 T=T+V( I Tl - 1 , IT2 —1,1 R-l )*A( IT2,IR)*X(IT2,IR) 

11 X ( I Tl, I R— 1 )=X( IT1.IR— 1)-T 

RETURN 
END 
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tIBFTC AE1 


FUNCTION AE( I, J.K) 

C 

C CALCULATION OF LEFT-HAND SIDE. 

C0MM0N/BL0CK1/NR.NT, RA! 17) t T A ( 25 ) , X ( 25 , 16 ) ,X1, AVS 

COMMON/ BL0CK2/V(25,25»15).A( 25, 15 )• FM( 25, 15) , BCMC ( 251 ,BCMR( 16) , 

1 BCMR1 (25) ,BCMA( 16 ) , VI ,FM1 , NRM1 ,NRM2, NRM3, NTM1 ,RBA ( 16) « DRA ( 16 ) , 

2 NTE.NTEMl.TBA(25).DTA(25)tSTBA(25)* ROW (25)«D(25.25) 

LOGICAL AVS 

C 

M = NT 

IF(I.EQ.l) GO TO 17 
IF(J.EQ.l) GO TO 6 
IF(J.EQ.M) GO TO 11 

C 

C REGION 1 (INTERIOR POINTS). 

IF(K.EQ.O) GO TO 5 
GO TO (1,2.3, 4), K 

L-COEFFICIENT. 

1 AE=(RBA( I-1)*(TBA( J)-TBA(J-1 ) ) >*(RBA< I - 1 ) * STBA ( J ) +RB A < 1-1) 

1 *STBA( J-l) )/(2.*DRA(I-l>) 

RETURN 

R-COEFF IC I ENT • 

2 AE= ( RB A ( I )*(TBA( J)-TBA( J-l) ) )*(RBA( I ) *STB A ( J ) +RBA ( I )*STBA( J-l) ) 
1 / ( 2 • *DRA ( I ) ) 

RETURN 

T-COEFFICIENT. 

3 AE=(RBA< I )-RBA( I-1))*(RBA< I ) *STBA < J-l ) +RBA ( I -1 ) *STBA ( J-l ) ) 

1 / ( 2. *R A ( I )*DTA< J-l) ) 

RETURN 

B-COEFF IC I ENT • 

4 AE= ( RB A ( I )— RBA( I — 1 ) ) * ( RB A( I ) *STBA ( J ) +RBA ( I -1 ) *STBA ( J ) ) 

1 / ( 2. *RA ( I > *DTA ( J ) ) 

RETURN 


D-COEFFIC IENT=L+R+T+B. 

5 AE=-(( RBA( I- 1 ) * ( TBA( J )-TBA( J-l ) ) ) * ( RBA (I -1 ) *STBA ( J ) +RBA ( 1-1 ) 

1 *STBA(J— 1 ) ) / ( 2. *DRA ( 1-1 )) 

2 +(RBA(I)*(TBA( J)-TBA( J-l) ) )*(RBA( I ) *STBA ( J ) +RB A ( I )*STBA( J-l) ) 

3 / ( 2 .*DRA ( I ) ) 

4 + ( RB A ( I ) -RBA (I— 1))*(RBA( I ) *STBA ( J-l ) +RBA ( 1-1 )* STBA ( J-l ) ) 

5 /(2.*RA( I)*DTA( J-l) ) 

6 + ( RB A ( I )— RBA ( I — 1 ) ) *(RBA( I ) *STBA ( J ) +RBA ( I - 1 ) # STB A ( J ) ) 

7 / ( 2 •* R A( I ) *DTA ( J ) )) 

' RETURN 


REGION 2 ( THETA=THETAO BOUNDARY LINE, J=l). 
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6 IF(K.EQ.O) GO TO 10 
GO TO (7,8,16,9) ,K 


L-COEFF IC IENT • 

7 AE=(RBA(I-1)*(TBA(1)-TA( 1) ) ) *( RBA( I-1)*STBA( 1) )/ (2.*DRA( 1-1 ) ) 
RETURN 

R-COEFF IC IENT • 

8 AE=(RBAtn*(T8A(l)-TA(l) ) ) * ( RBA ( I ) *STBA( 1) )/(2.*DRA(I ) ) 

RETURN 

B-COEFF ICIENT • 

9 AE=(RBA( I )-RBA(I-l))*(RBAf I) *STBAf 1 )+RBA{ 1-1 )*STBA( 1)1 
1 /(2.*RA( I ) *DTA ( 1 ) ) 

RETURN 

O-COEFFICIENT. 

10 AE=-((RBA(I-i)*(TBA( 1)-TA( 1) )>*(RBA( I -1 ) *STBA( 1 ) ) / ( 2 .*DRA ( 1-1) ) 

1 ♦(RBA(I)*(TBA(1)-TA(1) ))*(RBA( I )*STBA( 1) )/<2.*DRA( I ) ) 

2 + ( RB A ( I )-RBA( I-1))*(RBA{ I ) *STBA < 1 > +RBA ( I- 1 )*STBA ( 1 > ) 

3 /(2.*RA( I ) *DTA ( 1 ) ) ) 

RETURN 

REGION 3 (NEGATIVE AXIAL BOUNDARY, THETA=P I ) • 

11 IF(K.EQ.O) GO TO 15 
GO TO (12, 13, 14, 16), K 

L-COEFF IC IENT • 

12 AE=(RBA( I-1)*(TA(M)-TBA(M-1> ))*{RBA( I -1 ) *STBA { M-l ) )/(2.*DRA( 1-1) ) 
RETURN 

R— CDEFF IC IENT • 

13 AE= ( RB A ( I )»(TA(M) -TBA ( M- 1 ) ) )*(RBA( I ) *STBA( M-l ) ) / ( 2 .*DRA ( I ) ) 

RETURN 

T-C3EFF IC IENT • 

1A AE= ( RB A ( I )-RBA(I-l))»(RBA( I) *STBA( M-l ) +RBA ( 1-1 )*STBA ( M-l ) ) 

1 /(2.*RA( I)*DTA(M-1) ) 

RETURN 

D-CDEFF IC IENT • 

15 AE = { (RBA( 1-1 )*(TA( M)-TBA (M-l ) ) )*f RBA( 1-1 )*ST8A(M-1 ) ) 

1 /(2.*DRA( 1-1) ) 

2 *(RBA( I ) *( TA ( M )-TB A( M- 1) ) ) * ( RBA ( I ) *ST8 A ( M-l ) )/(2.*DRA(I) ) 

3 + (RBA ( I >-RBA(I-l ) ) * ( RBA(I) *STBA ( M-l ) +RBA ( I -1 ) *STBA ( M-l ) ) 

A /(2.*RA(I)*DTA(M-1) ) ) * 1-1. 1 

RETURN 

(NO B-COEFF IC IENT ) • 

16 AE=9. 

RETURN 

REGION A (CENTER OF SPHERE, 1=1). 

17 S= (RBA(1)*(TBA(1)— TA(1> ) ) * ( RBA( 1 ) *STBA (1) ) / ( 2 . *DRA( 1 ) ) 

*M1=M-1 

DO 18 J1-2,NTH1 

18 S=S+(RBA( 1)*(TBA( J1)-TBA(J1-1)))*(RBA( 1)*STBA( J1)+RBA( 1) 

1 *STBA( Jl-1) )/(2.*DRA(l) ) 

S=S+(RBA( 1)* (TA(M)-TBA(M-l)) ) * { R6 A ( 1 ) *STBA < M-l ) ) / ( 2 . *DRA ( 1 ) ) 

AE=-S 

RETURN 

END 
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$ IBFTC FE1 


FUNCTION FE ( I » J ) 

CALCULATION OF RIGHT-HAND SIDE. 

COMMON/ BLOCK 1/ NR , NT, RA(17).TA(25),X(25,16) ,X1,AVS 
C0MM0N/BL0CK2/V(25,25,15),A(25,15),FM(25,15>,BCMC<25),BCMR(16) 

1 BCMR1I25) f BCMA( 16 1 , V 1 , FM1 , NRM1 , NRM2 ,NRM3,NTM1 ,RBAI 16) ,DRA( 161 

2 NTE,NTEM1,TBA(25) ,DTA(25) , STBA ( 25 ) , ROW ( 25 > ,0C 25, 25 ) 

LOGICAL AVS 

M*NT 

IF(I.EQ.l) GO TO 3 
IF(J.EO.l) GO TO 1 
IF(J.EQ.NT) GO TO 2 

REGION 1 (INTERIOR POINTS!. 

FE = ( RB A ( I I-RBA (I— 1))*(R8A( I ) • ( TBA ( J ) -TBA ( J-l ) > +RBA( I-i ) 

1 *(TBA(J)— TBA(J-l) )) *(RBA( I)«(STBA( J)+STBA( J-l) )+RBA( 1-1) 

2 * ( STB A ( J ) +STBA ( J- 1 J ) )/8. 

RETURN 

REGION 2 ( THETA=THETAO BOUNDARY LINE. J = l). 

1 FE=(RBA(I )-RBA(I-l))*(RBA(n*(TBA(l)-TA(l) )+RBA(l-l) 

1 * ( TBA { 1 ) -TA ( 1 ) ) ) * ( ( RB A( n+RBAU-1) )*STBA(1) »/8. 

RETURN 

REGION 3 (NEGATIVE AXIAL BOUNDARY, THETA=PI ) . 

2 FE= ( RBA ( I I-RBA ( 1-1 > ) *(RBA( I I * ( TA ( M I -TBA( M-l ) )+RBA( 1-1) 

1 *(TA!MI— TBA(M-l) I )*{ ( RBA( I I +RBA( I-1I )*STBA(M-X) 1/8. 

RETURN 

REGION A (CENTER OF SPHERE, 1=1). 

3 S= ( TBA ( 1 1 -TA ( 1 1 ) *STB A ( 1 ) 

F Ml =M- 1 

00 A J 1 =2 , MM1 

A S=S+(TBA( J1)-TBA( Jl-1 » I * ( STB A ( J 1 1 +ST8A ( J 1-1 1) 
S=S+(TA(M)-TBA(M-1))*STBA(M-1> 

FE=RBA (11**3 *S/6. 

RETURN 

END 
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$ I BFTC FAC1 


SUBROUTINE F ACTOR ( A* ROW* OROER I 

SEE COMPUTATIONAL METHOOS OF LINEAR ALGEBRA 
BY FADEEVA (PAGE 86)* REFERENCE 2 • 

INTEGER ROW* ORDER 
DIMENSION A(25*25)#R0W(25) 

DOUBLE PRECISION T,ZERO 
REAL MAX 

DATA ZERO/O. ODO/ 

N=ORDER 

ASSIGN I TO K 1 
ASSIGN 9 TO K2 
DO 14 I *1 ♦ N 
IP1=I+1 
IM1=I-1 

GO TO Kl,(2,l) 

1 ASSIGN 2 TO K1 
GO TO 5 

CALCULATION OF LOWER TRIANGULAR MATRIX* 

2 DO 4 L = I * N 
T*ZERO 

DO 3 K= 1 * I Ml 

3 T=T+A(K,I )*A(L,K> 

4 A(L*I)=A(L*I)-T 

DETERMINE MAXIMUM ELEMENT IN COLUMN I. 

5 MAX=0.0 
DO 6 L = I * N 

IF(ABS(A(L,I)).LE.ABS(MAX))GO TO 6 
NI = L 

MAX = A( L * I ) 

6 CONTINUE 

IF»MAX.EQ.O.O)CALL ARERR (37H SUBR. FACTOR HAS A SINGULAR MATRIXES) 
ROW ( I ) - NI 

IFU.NE.NDGO TO 16 
IF< I.EQ.NJRETURN 
GO TO 8 

INTERCHANGE ROWS I AND NI. 

16 DO 7 K=1,N 
S=A ( I * K ) 

A ( I * K) - A ( N I * K ) 

7 A ( NI *K ) =S 

DIVIDE ELEMENTS OF ROW 1 
BY MAXIMUM ELEMENT A(l*l). 

8 GO TO K2» (11*9) 

9 DO 10 L=2 * N 

10 A(ltL) s A( 1 «L I/MAX 
ASSIGN 11 TO K2 
GO TO 14 

CALCULATION OF UPPER TRIANGULAR MATRIX. 

11 00 13 L = 1 P 1 1 N 
T^ZERO 

DO 12 K*= 1 * I M 1 

12 T=T+A( I * K ) *A(K*L ) 

13 ACI*L)=(A(I f L)-T)/MAX 

14 CONTINUE 
C 

RETURN 

END 
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$IBFTC INV1 


SUBROUTINE I NVERT ( A, ROW , GROER, D ) 

SEE COMPUTATIONAL METHODS OF LINEAR ALGEBRA 
BY FADEEVA (PAGE 87), REFERENCE 2* 

INTEGER ROW, ORDER 

DIMENSION A ( 25 , 25 ) ,ROW( 25)«D(25,25) 

DOUBLE PRECISION T,ZERO 

DATA ZERO/O. 0D0/ 

N*ORDER 

ASSIGN 7 TO Ki 
I*N 

1 IP1 = I + 1 
IM1=I-1 

USING LOWER TRIANGULAR MATRIX 
CALCULATE ELEMENTS OF L- I NVERSE. 

DO 4 M*l, I 
L- I Pl-M 
T*ZERD 

GO TO KI, (2,7) 

7 ASSIGN 2 TO KI 
GO TO 4 

2 LP1-L+1 

DO 3 K*LP1,N 

3 T=T+A(K,L)*D( I,K) 

4 D(I,L)-(FLOAT(L/I)-T)/A(L,L) 

IF(I.EQ.l)GO TO 8 

USING UPPER TRIANGULAR MATRIX 
CALCULATE ELEMENTS OF U-INVERSE. 

DO 6 M= 1 , I Ml 
L*I-M 
T*ZER3 
LP1=L*1 
DO 5 K*LP I , N 

5 T=T+A(L,K)*D(K,I) 

6 D(L,I)=-T 
I s I Ml 

GO TO 1 

INTERCHANGE COLUMNS I AND K. 

8 I*N-1 

9 K*ROW( I ) 

DO 10 L*1,N 
S*D(L,I) 

D( L , I ) *D( L,K ) 

10 D( L ,K) *S 

IF( I.EQ.DRETURN 
1*1-1 
GO TO 9 
END 
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APPENDIX C 


FORTRAN variable 
A(25, 15) 

AVS 


BCMA(16) 

BCMC(25) 

BCMR(16) 
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FORTRAN SYMBOLS 

Definition 

storage for diagonal matrices 

axial value switch (used in cases II and III as a logical 
quantity) 

boundary condition multiplier array on negative axis (0 = n) 

boundary condition multiplier array on surface of sphere 

boundary condition multiplier array on 6 = constant surface 
of segment 

boundary condition multiplier at central point 
matrix storage area 
Ar array 
A0 array 

F (right-hand side) multiplier array 
coefficient of center point 
looping index 

I - 1 
I + 1 

r indexes 
6 indexes 

miscellaneous switching integers 

looping indexes 

maximum value storage 
order of matrix 
matrix row index 
number of r mesh lines 


mu hi linn 


l l l 


i i i mill lll■lllllllll^lllll 


mu 


in mu mil ii in 


FORTRAN Variable 


Definition 


NRM1, NRM2, NRM3 

NR-1, NR-2, NR- 3 

NT 

number of 9 mesh lines for full sphere 

NTE 

number .of 9 mesh lines for spherical segment 

NTEM1 

NTE - 1 

NTM1 

NT - 1 

ORDER 

order of matrix 

RA(17) 

array of r mesh values 

RBA(16) 

r array 

ROW(25) 

matrix row indicators 

S 

summation storage 

STBA(25) 

sin 9 

T 

temporary storage 

TA(25) 

array of 9 mesh values 

TBA(25) 

9 array 

V(25, 25, 15) 

25 by 25 matrix storage area; stores up to 15 inverted 
matrices 

VI 

inverse of center point 

X(25, 16) 

interior mesh values 

XI 

central mesh value 

Y(25) 

double precision temporary storage array 

ZERO 

double precision initializer 



APPENDIX D 


FLOW CHARTS 

In the following flow charts (figs. 12 to 21) all diamond shapes refer to decision mak- 
ing points in the subroutines. All other operations are indicated with rectangular blocks. 
Numbers appearing on the upper left of any block refer to the corresponding statement 
numbers in the program. All FORTRAN symbols used are defined in appendix C. 



Figure 12. - Macro flow chart showing interrelationship of various subroutines 
used in program. l !, Uses functions AE and FE. ) 
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Figure 14. - Subroutine SPHER1. 



























Figure 15. - Subroutine SPHER2. 



Figure 16. - Subroutine SEG1. 












Figure 17. - Subroutine SEG 2 . 





Figure 18. - Function AE(I, J f K) - calculation of LHS coefficients. (Formulas for 
quantities L, R, T, B, and D are given in the main text.) 


Interior points 


Figure 19. - Function FE(I, J) - calculation of RHS coefficients. (See text for formula for F.) 
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TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech BriefSj 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



